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We discuss the implementation ol optical quantum networks where the interface between sta- 
tionary and photonic qubits is realized by optomechanical transducers [K. Stannigel et al, PRL 
105, 220501 (2010)]. This approach does not rely on the optical properties of the qubit and thereby 
enables optical quantum communication applications for a wide range of solid-state spin- and charge- 
based systems. We present an effective description of such networks for many qubits and give a 
derivation of a state transfer protocol for long-distance quantum communication. We also describe 
how to mediate local on-chip interactions by means of the optomechanical transducers that can 
be used for entangling gates. We finally discuss experimental systems for the realization of our 
proposal. 
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I. INTRODUCTION 

The distribution of quantum information between in- 
dividual nodes of larger quantum networks is a key re- 
quirement for many quantum information applications 
PP, in particular for long-distance quantum communi- 
cation and quantum key distribution protocols [2]. In 
contrast to classical networks, quantum communication 
channels must allow the distribution of entanglement, 
which is the essential resource that can be harnessed 
by further local processing within the nodes. In build- 
ing high-fidelity quantum channels over long distances, 
qubits encoded in propagating photons play a unique 
role, since they provide the only way to transmit quan- 
tum information over kilometer distances [3J 0]. How- 
ever, also on a smaller scale rapid progress is made in 
the design of nano-fabricated photonic circuits for classi- 
cal and quantum applications [5] and compared to elec- 
tric circuits [7] or diverse qubit shuttling techniques 
[3 |9], such photonic circuits could provide a fast and 
robust alternative for 'on-chip' distribution of entangle- 
ment. Therefore, the development of coherent interfaces 
between stationary and so-called 'flying' photonic qubits 
is one of the essential steps in realizing quantum net- 
works and is commonly included in the list of key crite- 
ria for general purpose quantum information processing 
platforms [TU] . 

The fundamental operation to be carried out in a quan- 
tum network is the transfer of an arbitrary quantum state 
between two nodes according to |V>)i|0)2 — > \0}i\ip}<2, 
which together with local operations allows the genera- 
tion of inter- node entanglement. Ideas for a physical im- 
plementation of this basic building block have first been 
developed in the context of atomic cavity QED. Here, 
qubits are stored in spin or hypcrfine states of trapped 
atoms or ions and can be selectively excited to other elec- 
tronic levels where they couple strongly to light. As was 
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FIG. 1. (Color online) (a) Elementary quantum network for 
the realization of state transfer protocols between two nodes 
via the exchange of photons propagating along an optical 
fiber. Tunable qubit decay rates R allow for a controlled 
emission of the photon wave-packet and a perfect reabsorp- 
tion at the second node, (b) Schematic setup for implement- 
ing effective qubit-light interfaces based on optomechanical 
transducers (OMTs): The mechanical resonator mediates a 
coupling between a solid-state qubit and the driven optical 
cavity mode. As a result, qubit excitations decay with an ef- 
fective rate T into the fiber via the cavity output (see text for 
details). 



shown by Cirac et al. [IT] , this level of control can be used 
to design a deterministic state transfer protocol based on 
the tailored emission of a photon, which is subsequently 
reabsorbed at a second node along an optical fiber (see 
Pig-Ba)). 

In parallel to the developments in the field of atomic 
qubits, substantial progress has been made on nano- 
fabricated solid-state qubits, including for example 
impurity-spins [12] in various host-materials (such as sil- 
icon P2] or diamond [HI [T5]). quantum dots [THl [T7] 
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or superconducting devices [TSH2T] . In view of the re- 
markable level of coherence and control that has already 
been achieved in such systems, the challenge is now to 
identify a suitable interconnection with optical quantum 
channels also for this larger class of qubits. This is, how- 
ever, often hindered by lack of coherent optical transi- 
tions or the incompatibility with light. Early suggestions 
to overcome this problem have been made in the context 
of hybrid systems, where, e.g., superconducting devices 
are coupled to close-by atomic [22H24] , molecular [25] [26] 
or spin systems [2TH21?] . An optical interface could then 
in principle be realized by using, for example, the atoms 
as a mediator. Recently, an alternative promising route 
has become available which is fully solid-state based and 
makes use of the quantized motion of macroscopic me- 
chanical resonators. In a cavity-optomechanical setting, 
these resonators can coherently interact with light [30] 
and their control near the single-quantum level has been 
demonstrated [3TH51] . On the other hand, they can also 
interact with solid-state qubits via magnetic [3"5H3"T] or 
electric fields [55HS] , and these abilities make them nat- 
ural candidates for mediating indirect qubit-photon in- 
teractions, thereby realizing the desired qubit-light in- 
terface. 

In a recent work [32] we proposed the use of such 
optomechanical transducers (OMTs) for the implemen- 
tation of optical communication protocols between two 
solid-state qubits, while related ideas have also been dis- 
cussed in the context of traveling- wave phonon to photon 
converters [33] ■ The basic idea of an OMT is illustrated 
in Fig.fiTb), where a spin- or charge-based qubit is cou- 
pled to the motion of a mechanical resonator via mag- 
netic field gradients or electrostatic interactions. The 
resonator in turn interacts with the field of an optical cav- 
ity mode via radiation pressure or optical gradient forces 
as is currently experimentally explored with various dif- 
ferent optomechanical (OM) settings [3JH31 E2H3S] ■ As 
a result, this configuration induces effective interactions 
between the qubit and photons which do not rely on the 
optical properties of the qubit system. Therefore, this 
approach in principle enables quantum communication 
applications for a large range of solid-state qubits. 

The purpose of the present work is two-fold. In 
the first and main part of the paper we study long- 
distance optomechanical quantum networks as intro- 
duced in Ref. [32]. Here, we present a detailed analysis 
of the effective qubit-fiber interactions and derive a cas- 
caded master equation which describes the effective mul- 
tiqubit dynamics in such a network. In particular, we 
describe the implementation of tunable qubit-fiber cou- 
plings for quantum state transfer protocols and discuss 
the role of thermal noise and Stokes-scattering processes. 
These are the two main noise sources which are typi- 
cally absent in analogous atomic settings. In the second 
part of the paper we then investigate the case of a single 
mode optical quantum channel relevant for optical on- 
chip communication schemes. In this setting, the qubits 
are effectively connected by a finite number of mechani- 



cal and optical modes and we analyze the performance of 
entangling operations mediated by this type of network. 

The remainder of this paper is structured as follows: 
In Section |n] we present a detailed description of a single 
OMT as depicted in Fig.jljb) and discuss the resulting 
effective qubit-fiber interaction together with the rele- 
vant noise processes in this system. In Sec. |III| we extend 
the model to an arbitrary number of nodes and derive 
a cascaded master equation for this multiqubit network. 
Further, we give a protocol for state-transfer between 
two nodes and assess its performance. In Sec. IV we 
then study an on-chip setting where the individual nodes 
are coupled by optical resonators with a discrete spec- 
trum and analyze the prospects for performing entangling 
gates between the qubits. Finally, we discuss potential 
experimental realizations of the proposed OMT in SecfV] 
and close with concluding remarks in Sec.|VI| 



II. AN OPTOMECHANICAL INTERFACE 
BETWEEN STATIONARY QUBITS AND LIGHT 

The aim of this work is to develop a theoretical 
model for the description of optical quantum networks 
where the interconversion between stationary and pho- 
tonic qubits is realized by an OMT at each node. As a 
first step we will investigate in this section the dynamics 
of a single node of such a network and derive a simplified 
model for the effective qubit-light interactions mediated 
by the OM device. 



A. Model 

We consider a single node of an optical network as 
schematically shown in Fig.[ljb), where the coupling be- 
tween a qubit and an optical fiber is mediated by an 
OMT. The total Hamiltonian for this system is 



H = H n 



>dc 



lib 



-fib 



H, 



(1) 



where iJ no dc describes the coherent dynamics of the qubit 
and the OM device, the free evolution of the optical 
fiber modes, and -ff C av-fib the coupling between cavity and 
fiber. The last term in Eq. 0, H cnV7 summarizes all 
additional interactions with the environment leading to 
decoherence as specified below. 

For a general discussion we assume complete local con- 
trol over the qubit, which is taken to be encoded in two 
states |0) and |1) of an isolated quantum system, and 
which interacts with the quantized motion of a mechan- 
ical resonator, e.g., via magnetic field gradients |36) or 
electrostatic interactions [58H1T] . The mechanical res- 
onator in turn is coupled to the field of an optical cavity 
mode via radiation pressure or optical gradient forces. 
In terms of the usual Pauli operators a z = |1)(1| — |0)(0| 
and a~ — |0)(1| for the qubit, and the bosonic operator 
b for the mechanical mode, the Hamiltonian H„ n Af, can 
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be written in the form (ft = 1) 

-ffnodc = + ^{a-tf + a+b) + H om . (2) 

Here, uj q is the qubit's energy splitting, A the strength of 
the qubit-resonator coupling and H om is the Hamiltonian 
for the coupled OM system. We point out that depending 
on the physical implementation, the Jaynes-Cummings 
type interaction assumed in Eq. |2]) may emerge only as 
an effective description of an underlying driven two- or 
multilevel system. In this case, also Lo q is an effective 
parameter which can differ from the bare frequency scales 
in the system. (This will be discussed in more detail 
in Sec.|V] where we consider the implementation of the 
Hamiltonian pi) for the specific examples of spin and 
charge qubits.) 

The OM system consisting of a mechanical resonator 
coupled to a single optical cavity mode is described by 
the Hamiltonian 

H OI[L — to r b^b + w c c^c + <7o c c (b + tf) (3) 
+ iE(t)e- WLt c) -i£*{t)e iujLt c, 

where u> r is the mechanical vibration frequency and c is 
the bosonic annihilation operator for the optical mode 
of frequency uj c . The quantity go = ao^r is the single 
photon OM coupling and corresponds to the optical fre- 
quency shift per zero point motion a^. In the last line 
of Eq. (§ we have included an additional external laser 
field of frequency Wi, which coherently excites the cav- 
ity field. As detailed below, this driving field will allow 
us to enhance and also to control the OM coupling by 
adjusting the slowly varying driving strength £ (t). We 
finally stress that the optical cavity might support ad- 
ditional degenerate modes which we ignored in writing 
Eq. As will be discussed below, this is valid as long 
as the OM interaction does not cause scattering between 
the degenerate modes. 

B. Quantum Langevin equations 

The optical cavity of the OM system is coupled to the 
modes of an optical fiber, which in the multinode set- 
ting considered below will serve as our quantum commu- 
nication channel. The fiber supports many continua of 
modes described by bosonic operators a au] , normalized to 
[a aa) , dg-i^i] = S ayCr 'S(bj — a/), where a labels, e.g., propa- 
gation direction or polarization. The free fiber Hamilto- 
nian is then simply given by Hgb = J2a I w aJ. w a CTW . 
We assume that the OM system is located at position 
z = along the fiber to which it is coupled according to 

#cav-fib = j ^ (V K <xM ca L - h - C -) I 

where the coupling constants y K a (oj) are defined in this 
way for later convenience. For definiteness, we consider 



the case where the cavity mode is side-coupled to a fiber 
with left- and right-moving modes of a single relevant 
polarization (a — R,L). In this case, the cavity-fiber 
coupling is naively estimated to be proportional to the 
overlap integral 



VWM~ / dz'h(z')e^'/ c , (4) 

J -Az 

where h(z) is the mode function of the evanescent cav- 
ity field over the relevant interaction region of length 
2Az (see, e.g., Refs. [13 US] for a detailed analysis of 
the coupling to a whispering-gallery mode cavity). For 
Az > A c , where A c is the vacuum cavity wavelength, 
we can roughly distinguish between two cases. For a 
standing- wave cavity mode h(z) ~ cos(w c z/c) and pho- 
tons from the cavity will be emitted into both direc- 
tions (k,r ~ k-l), while for a running- wave mode such 
as preferred emission into a specific di- 

rection can be achieved (in this case kr ^> «&). Although 
in principle both configurations could be relevant for dif- 
ferent quantum communication applications, we will for 
concreteness focus in the following on the case where c 
describes a circulating running-wave cavity mode that 
couples to the fiber according to 

tfcav-fib - (cf R { Z = 0) - cUr{z = 0)) , (5) 

where fn(z) = -|= J °° duj aji U) e lu}Z / c is the right- 
propagating fiber-field. The coupling of the cavity to this 
continuum of modes will lead to an irreversible decay of 
the cavity field at the rate k/ = kr(uj c ). 

To proceed, we eliminate the time-dependence in 
Eq. (J3j) by transforming all photonic operators to a frame 
rotating at the laser frequency. This amounts to replac- 
ing H om by 

H om = u r tfb + A c c f c + .g c f c(6 + 6 f ) (6) 
+ i(£(t) C t-£*(t)c) , 

and Hfi h by 5 fib = J2c I duj ^ a lu a <?u , where A c = 
uj c — ujl and A w = oj — lo^ are the detunings of the cavity 
and fiber modes from the laser frequency, respectively. 
We then use a standard Born-Markov approximation to 
eliminate the fiber modes [49] and describe the resulting 
dissipative dynamics in terms of a quantum Langevin 
equation (QLE) for the cavity field, 

CKi-i[c, -ffnodc] - KC- y/2Kff in (t) - \/2K^/o(*) . (7) 

Here, we have introduced a total decay rate k = Kf + 
kq, where kq accounts for additional intrinsic losses of 
the optical cavity mode and fo(t) is the associated noise 
operator. We have further defined operators fi n (t) :— 
/fl(0~,t) and / ou t(t) := /i?,(0 + ,t) which obey the input- 
output relation 

/out(0»/in(*) + >/2^c(t), (8) 
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and we specify vacuum white noise statistics for the fiber 
input, i.e., [f in (t)jl(t')} = 5(t-t>) and (f in (t)fl(t')) = 
S(t — t'), as well as for fo(t). In addition, the mechanical 
resonance generally has a finite intrinsic width and is con- 
nected via its support to a thermal bath of temperature 
T. This is commonly modeled by the QLE 



-i[b,H a 



7r> 



(9) 



where 7 m /2 is the decay rate of the mechanical ampli- 
tude. In contrast to the optical fields, the mechanical 
noise is characterized by a thermal white noise opera- 
tor £(i) with statistics = ( N m + l)S(t - *') 
and [£(i),C f (i')] = S(t - f). Here, N m is given by the 
Bose occupation number for a mode of frequency uj r . In 
the high-temperature case, where fceT ^> Htu r and hence 
N m 3> 1, the relevant thermal decoherence rate will turn 
out to be r y m N m w where Q — (jj r /j m is the quality 
factor of the mechanical resonance and typically Q 1. 
Below, we will often take the limit 7 m — > and whenever 
doing so it is understood that we keep the decoherence 
rate j m N m constant. 



C. Linearized equations of motion 

In typical experiments the OM coupling go is too small 
to allow coherent interactions between the mechanical 
system and individual cavity photons. Therefore, to 
achieve an appreciable and also tunable coupling we con- 
sider the case of a strongly driven OM system where the 
effective coupling between mechanics and light is ampli- 
fied by the coherent field amplitude inside the cavity. 

Starting from the QLEs ([7]),([9|, we perform a unitary 
transformation c — > c + a and b — > b + j3 such that the 
c- numbers a, (3 describe the classical mean values of the 
modes and the new operators c and b represent quantum 
fluctuations around them. We require any classical (that 
is c-number) contributions to the transformed QLEs to 
vanish, which yields the following system for the classical 
response: 

a = -(iA c + K)a - ig {fi + f3*)a + £(t) , (10a) 

$ = -(W r + 7m/2)0 " ig M 2 . (10b) 

For £{t) = const and j m — > the steady state solution is 
(3 = —go\a\ 2 /uj r with a determined from 



£ 



a 



iA c + k — 2i g^ \a\ ju) r 



(11) 



This relation is still approximately valid if £(t) is slowly 
varying compared to the characteristic response time. 
This means, e.g., £(t)/£(t) -C K, A c) oj t for the case 
\a\ go -C A c , uj t which is of relevance below. Within these 
limits, a desired temporal profile a(t) can be directly re- 
lated to the applied laser power and phase by means of 
Eq. (111. For more rapid variations the system (10 1 has 



In the strong driving regime where \a(t)\ 2 > 1 the OM 
coupling of Eq. Q can be linearized: After going to the 
displaced operators in the QLEs ([7|,([9| and eliminating 
the classical forces, we drop terms of order go relative 
to those of order go\a\ and go \a\ 2 . The result of this 
procedure is equivalent to replacing H om by the linearized 
OM Hamiltonian 



u) r tfb + A c (t)Jc+(G(t)J 



G*(t)c)(b + tf). 

(12) 



Here, we have introduced the laser-enhanced OM cou- 
pling G{t) = gou{t), which describes interconversion 
between phonons and photons, as well as the renor- 
malized cavity detuning A c (t) = A c + g ((3 + 13*) — 
A c — 2 |G(t)| 2 juj r . Although the latter shift is small 
in the regime of interest identified below, one has to 
take it into account when varying G(t). We will drop 
the tilde in what follows and also note that we have ig- 
nored a contribution X/2(/3(t)a + +h.c.) to H DO d c , which 
can be compensated by local fields acting on the qubit. 
Finally, if the cavity supports many degenerate modes 
c M that couple to the mechanical resonator according to 
9o Yin c li c v(b + tnGn on ly t ne OM coupling to the 
driven mode is enhanced. The others may then be ne- 
glected due to go \G\, which justifies our single-mode 
treatment. However, this argument breaks down if there 
is phonon-assisted scattering between the modes, i.e., for 
an interaction of the form go v c^c u {b + 

In summary, the OM system can be described by an 
enhanced, linear coupling \G\ ^> go in ffjj between the 
resonator and an undriven cavity mode. The qubit dy- 
namics is determined by the Heisenberg equations of mo- 
tion 



_ A , 

a = —iuj q a +i — a z b, 

a z = -i\(a + b-tf<7-), 



(13a) 
(13b) 



and the linearized OM dynamics (QLEs Q, ^ with 
H om — > i?om) can t" e written compactly by introducing 
the vectors v = (6, c, b\ c^) T and S = (tr _ , 0, — cr + , 0) 



T. 



v(t) = -Mv(t) - i~S(t) - N(t) . (14) 

Here, the drift matrix M describes the response of the 
OM system and is explicitly given by 



M 



1 2 

G 


-CG* 



G* 
A r — in 







-CG* -u r - 



CG 


-G 







-G* 



-A, 



(15) 



to be integrated exactly to capture all retardation effects. 



where C = 1 gives the full linearized OM coupling 
and C — corresponds to the often-applied rotating 
wave approximation (RWA), which is valid for K,\G\ <C 
A c uj r and renders the coherent part of the dy- 
namics excitation-number- or energy-conserving. The 
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noise sources driving the OM system are summarized in 
N(i) = v^I(i) + R(t), where 1(f) = (0, / in , 0, // n ) T 
represents the fiber input-field and the intrinsic noise is 
contained in R(t) = (^7^£, V^o/o ,VTW^i V / 2ko/o) T - 
The statistics of these noise vectors follow directly from 
the definitions in the previous subsection and can also 
be found in App.|Bj We finally note that the system de- 
scribed by M generally exhibits self-oscillations for blue 
detuning A c < as well as bi-stability for red detuning 
A c > and strong coupling \G\ 2 > (k 2 + A 2 )w r /4A C 
(see, e.g., Ref. [50]). However, in this work we will only 
consider parameters where these phenomena do not oc- 
cur. 



D. Effective qubit-light interface 



The dynamics of a single node as described by Eqs. ( 13 ) 



and ( 14 ) allows for a coherent conversion of a qubit ex- 



citation into a photon propagating in the optical fiber, 
as is indicated in Fig.[2ja): First, the excitation is trans- 
ferred to the mechanical resonator at a rate ~ A, from 
where it is then up-converted by the excitation-number 
conserving terms of the OM interaction into a cavity pho- 
ton with a rate ~ G, and subsequently emitted into the 
fiber at a rate ~ Kf. At the same time, this transfer is 
affected by decoherence in form of mechanical dissipa- 
tion and photon loss, thereby degrading the fidelity of 
the interface under realistic conditions. Our goal is now 
to derive a simplified model for the interface which cap- 
tures all of these aspects in an effective description for 
the qubit-fiber coupling only. 

To proceed, let us first look at the dynamics of the 
linearized OMT, which is fully characterized by the drift 
matrix M defined in Eq. (15 1 (see also Ref. [5T] for an 
extensive discussion). It has two independent eigenval- 
ues which we identify with iuj± + j± , where w± > are 
the frequencies and j± the decay rates of the associated 
OM normal modes (the remaining eigenvalues are simply 
complex conjugates). We take G to be real and positive 
for simplicity and focus on the case A c « u> r relevant for 
resonant excitation transfer, where we can distinguish 
between two regimes: For weak OM coupling, G < k/2, 
the resonator and cavity modes retain their respective 
characters and we have, e.g., cj_ w w r , uj + A c . How- 
ever, due to a small admixing between the modes the 
optical cavity provides an additional decay channel for 
the mechanical system and can lead to significant mod- 
ifications of the mechanical damping rate. In the limit 
G <C k it is given by 7_ = -f, n /2 + G 2 / n (see left panel of 
Fig.|2jb). In the opposite regime of strong OM coupling, 
G > k/2, resonator and cavity hybridize and the OM 
system exhibits a normal mode splitting (uj± ss uj r ± G) . 
This hybridization also implies that both modes decay 
with the same rate y± a k/2 (right panel of Fig.[2jb)). 
In more general terms we can introduce an OM damp- 
ing rate 7 op , which for arbitrary parameters is defined 
a s 7op = m i n {7+>7-} an d reproduces 7 op ~ G 2 /k and 



7 op a k/2 in the two respective limits. As 7 op is the 
rate at which mechanical excitations are converted into 
traveling photons, it sets a maximal transfer rate for the 
OMT. We expect 7 op 3> 7 m /2 for typical parameters, 
which justifies taking limit 7,„ — > 0. 

We now consider the experimentally relevant regime 
where the qubit-resonator coupling A is small com- 
pared either to 7 op or to the detuning from the near- 
est normal mode, that is, the regime where A <C 
max(7 op , min(|o;± — Wq|)). In this limit, the dynam- 
ics of the OM system is fast compared to the qubit- 
resonator coupling and can hence be adiabatically elim- 
inated, yielding the desired effective description of the 
OMT. To this end, we formally integrate Eq. (14) to ob- 
tain 



v(t) = v frco (i) - i\ [ dse-^*-)S(s) . 



(16) 



where we have neglected transients and Vf ree (f) = 
- f ^ dse~ M (*~ s )N(s) is the OM steady state solution 
in the absence of the qubit (see Eq. ( 14 ) for A —¥ 0). We 
may now in a Born-Markov approximation and to low- 
est order in A replace the qubit operators by their free 
evolution, i.e., cr~(s) = o~ (t)e lLUqt - t ~ s ^ and re-insert the 
result into Eq. ( 13 1. By keeping only the resonant terms, 



we obtain effective QLEs for the qubit operators, 



i(u} q + A ) + - 



a- - VTa z F in (t) , (17a) 



-r(l + o z ) + VAT (a+F in (t) + FiL(t)a-) , 



(17b) 



with noise operators i*In(t) = —iyj A 2 /4r 6f rcc (f). From 
Eq. pj) we further obtain an effective input-output rela- 
tion for the fiber field, 



fo 



■qTe^a-it). (18) 



In Eqs. ( 17 )-( 18 1 we have introduced the effective decay 



rate T and frequency shift Ao by the relations 
r=yRe{A u «)}, A = jIm{A n ( w ,)} , (19) 

where A(lS) = (M — itut)^ 1 is the OM response matrix. 
Further, rj = Kfjn is the branching ratio resulting from 
intrinsic cavity decay and ef> = a,rg{iA2i(ujq)}. Note that 



in writing Eq. ( 18 ) we have neglected small corrections 
related to finite j m , as well as non-RWA corrections and 
counter- rotating terms oc a + {t). 



The effective equations of motion (17 1 together with 



the input-output relation ( 18 1 arc familiar from atom- 
light interactions [49] and describe a two level system 
which is directly coupled to an optical fiber with an ef- 
fective decay rate Y. The OM system which mediates this 
interaction has disappeared from the dynamics, but still 
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FIG. 2. (Color online) (a) Level scheme illustrating the decay of a qubit excitation into the fiber (only processes corresponding 
to £ = are show). The dashed box represents the OMT and wavy arrows indicate mechanical diffusion and intrinsic cavity loss, 
(b) Level structure of the OMT as seen by the qubit. Wavy arrows indicate decay into the optical fiber at 7 op w min{G 2 /fv, k/2} 
in the regimes of weak (left) and strong (right) OM coupling, (c) Effective qubit-decay T as a function of the linear OM coupling 
G and the qubit frequency ui q as given by Eq. ( 19 1. We have used k = 0.05cj r and y m = 0, and A c = ui r is satisfied at G = 3k/2. 
(d) r as a function of cavity detuning and qubit frequency, where G = 3k/2 at A c = ui r . The dashed lines in (c),(d) indicate 
possible paths for tuning Y between ~ and ~ A 2 /2k (see text). 



determines T, Ao, <j>, which can thus be controlled by ap- 
propriately adjusting OM parameters. This can be seen 
in more detail from the expression for the decay rate, 



while the qubit frequency is fixed at ui q 



3k/2 



A 2 G 2 k/2 



(G 2 + (A c — u q )(uj q — uj r )) 2 + n 2 (uj q — uj r y 



(20) 



where we have assumed £ = and 7 m — > 0. A plot of the 
full expression from Eq. ( 19 ) is given in Fig.[2|c) for the 
case uj t rs A c and one clearly observes the underlying 
mode structure of the OM system as discussed above. 
The decay Y is suppressed if the qubit is detuned from 
the OM eigenmodes, while under resonance conditions, 
w, ~ w±, it scales as Y ~ A 2 /7 op , and we obtain Y w 
A 2 k/2G 2 and Y ps X 2 /2k in the weak and strong coupling 
regime, respectively. 

As we will discuss in more detail in Sec. |III C"] below, it 
is essential for the implementation of deterministic quan- 
tum state transfer protocols within quantum networks to 
have control over the shape of the emitted photon, which 
can be accomplished by considering a time-dependent de- 
cay rate Y — > Y(t). To achieve this, we may make use of 
the fact that in principle any of the OM parameters en- 
tering r could be varied in an experiment, that is, we 
consider Y{t) = Y(G(t), A c (t), iv q (t), ...). The effective 
Markovian description of the qubit dynamics and the al- 
gebraic relation in Eq. pl)| ) remain valid as long as the 
OM parameters are tuned slowly relative to the time- 
scale of the effective resonator decay l/7 op , such that 
G/G -C 7o P , etc. In this sense we may adiabatically 
adjust the effective qubit decay rate and hence acquire 
control over the generated signal photon. 

Examples. As a specific example we consider tuning 
the OM coupling along a path as indicated by the dashed 
line in Fig. gc). Here, G(t) varies between and ~ 3k/ 2 



and ui r ~ A c . This corresponds to tuning the lower nor- 
mal mode in and out of resonance with the qubit, such 
that r can be varied within a certain range [r m ; n ,r max ]. 
For G — > the qubit is far detuned by S = \u) g — w r \ 
from the resonator mode to which it couples and there- 
fore, the effective decay rate is suppressed according to 
T - A 2 7 op /2r5 2 « A 2 G 2 /2 K r5 2 , which realizes T min 0. 
A residual decay is still provided by mechanical damp- 
ing F m i n w A 2 7 m /4rJ 2 , which, however, is negligible for 
the parameters considered below. By increasing G up to 
~ 3k/2 we reach the opposite limit of strong OM cou- 
pling, where the lower OM normal mode is in resonance 
with the qubit and a maximal decay of r max = A 2 /2k is 
achieved. Within these limits, a desired time-dependence 
of the effective decay Y(t) can thus directly be translated 
into a profile Git) via Eq. (19). 



Similarly, we can vary the cavity frequency along the 



dashed line in Fig.[27d), where we chose 



-3re/2 



and it is assumed that the driving laser is at constant 
power. This means that G varies according to Eq. ( |11[ ), 
which has been included in the plot. We choose G = 3k/2 
at A c = oj r , such that r max for this path is the same 
as for the previous one. Starting from this point, the 
effective decay can be turned off by increasing the cavity 
detuning as far as needed, with a residual decay due to 
the resonator as before. 



E. Noise and imperfections 

As already mentioned above, the OMT does not only 
mediate the desired coherent interactions between the 
qubit and the optical channel, but under realistic con- 



ditions also adds noise. In the effective QLEs(17) this 



7 



noise appears first of all in form of the operator Fi n (t), 
which generally represents an effective non- vacuum input 
that heats the qubit. Second, noise processes also affect 
the out-field f ut(t) given in Eq. (18 1 where on the one 
hand the term Cf rcc (t) represents excess photons which 
are emitted into the fiber independently of the qubit state 
and on the other hand r\ < 1 accounts for loss of qubit 
excitations through the cavity's intrinsic decay channels. 

Let us first look at the noise operator i 7 i n (t) which ap- 
pears in the QLEs (JTYJ) for the qubit. From its definition 
in terms of the OM drift matrix M we see that two-time 
correlations of Fi n (t) decay with 7 op ^> T. Therefore, on 
the scale of the effective qubit dynamics, it can be con- 
sidered as nearly white noise and is characterized by an 
effective occupation number 

A = 2Rc / dt{fl(t)F in (Q))e- iw < t . (21) 
Jo 

This expression can be readily evaluated in Fourier space 
as shown in App.[B]and for j m — > we obtain 

Nq w 7miVm K 2 + (Ac - Uq? | K 2 + {A C -U q ) 2 



2 k 



G' 2 



4A c w 9 



The first contribution (given for £ = 0) is due to the 
thermal noise of the resonator's environment and roughly 
scales as "f m N m /^f op . The second contribution to Ao de- 
scribes heating due to non-energy-conserving terms in the 
OM coupling, ~ Gc'b\ and accounts for photons scat- 
tered from the strong driving beam into the cavity while 
simultaneously creating a mechanical excitation. We see 
that for a qubit on resonance with an OM normal mode 
(bj q « A c ± G) the conditions for small thermal noise 
are 7 m iV m <C 7op, while to reduce the Stokes scattering 
events we need k, G <C A c ,uj r . The latter conditions de- 
scribe the so-called resolved side-band regime and there- 
fore, we conclude that the conditions for a low noise OMT 
are equivalent to the requirements for OM ground state 
cooling [3"T] . 

The second type of noise is the contamination of the 
cavity output with noise photons originating from up- 
converted thermal noise or Stokes scattering events, as 
described by the term cx Cf lee (t) in Eq. (18). To ana- 



lyze the effect of these excess photons let us in a first 
step imagine that at time t = the qubit is prepared 
in the excited state and that we record the number of 
photons emitted into the fiber during a time-interval 
T ~ 1/r. The observed photon number is given by 

N(T) — f Q dt (/o U t/out) (t) and for simplicity, we assume 
that the qubit operator evolves according to cr~(t) = 



a (0), while being uncorrelated with Cf roc (i). 



We then obtain N(T) = 1-e 
ber of excess photons given by 



-iV ex (T) with the num- 



N m (T) = 2 Kf / dt(cl cc c itcc )(t), 



(23) 



which simply yields N CX (T) — 2k fT (c£ ree Cf ree )(0) for 
time- independent OM parameters. For T = kq = 0, 



7m — > 0, and u r = A c this quantity can be estimated [] 
to be A^r- 1 ) » 7m iV m /r + {n/T){G 2 /u*), and for a 
qubit resonant with an OM mode we may use T ~ A 2 /7 op 
to obtain the scaling 



A^r- 1 ) 



7 2 
I op 

A 2 " 



7> 



7op 



_K_G_ 

7o P w 2 



(24) 



While the terms in parentheses are small under the low 
noise conditions identified above, the prefactor is much 
larger than one due to the condition 7 op 3> A underpin- 
ning the adiabatic elimination. The number of excess 
photons may thus exceed one and make the signal pho- 
ton difficult to detect in general. However, this estimate 
did not take into account the fact that the qubit emits 
in a bandwidth V around uj q , while the OM noise floor 
is spread over a bandwidth 7 op . Therefore, we expect 
that the excess photons can be substantially reduced to 
A r ex (r _1 ) r/7 op by appropriate filtering, which cancels 
the large factor 7o P /A 2 appearing in Eq. (24). Indeed, 



we find below that the OMT at a second node provides 
exactly this filtering. 



III. QUANTUM NETWORKS 

With the concept of the OMT laid out in the previous 
section we are now ready to consider a quantum network, 
where the OMTs in the various nodes serve to link the 
qubits to a common optical fiber, as depicted in Fig.[T] 
To obtain an effective description of the multiqubit dy- 
namics, we proceed as before by linearizing the OM cou- 
plings and eliminating the coupled OMTs. We will then 
discuss how to realize a state-transfer within a two-node 
network. The Hamiltonian for the complete network is a 
straight-forward generalization of Eq. and reads: 



H 



node 



H 



cav-fib 



HI 



lib ■ 



(25) 



where H^ odc is given by Eq. ^ with appropriate node 
indices added and H l cnv summarize decoherence effects. 
We assume that the output of a given cavity is routed 
to the next, which is, e.g., realized by side-coupling all 
cavities to the right-moving field in the fiber as described 
by 



HI. 



-fib 



24(^(^-4/^)). (26) 



Here, /fl(zi) is the right-moving field operator defined 
after Eq. ^ and the nodes are located at positions < 



1 The steady state cavity occupation can be expressed 
as (ctc) free = f f^rCja M, with the matrix C J1 (w) 
defined in App.[B] Neglecting non-RWA corrections 
to the thermal contribution, we obtain (c^c)f rcc ~ 
I ^ bmN m \A 21 {u))\ 2 + 2K\A 2i (u)\ 2 ) which we evaluate 
to leading order in G 2 /u 2 . |52| . 
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Zi+\ along the fiber, with c$ denoting the driven cavity 
mode at node i. 

In a first step, we eliminate the modes of the optical 
fiber in a procedure similar to the one in Sec. |II B| For 
each node, we obtain an input-output relation of the form 
^ and in addition, we find a cascaded coupling where 
the output of a given cavity drives the subsequent one 



Ci(t) W i[Hl ode ,Ci] - KfCi(t) - V2Kff illt i(t) . 



where v* 



(hi, Ci, b\, c\) T and M l is defined as 



fin,i(t) ~ fout,i-l(t 



-i e 



Here, H l OYn in H^ odc is replaced by the analogue of Eq. 
Tij = (zi — Zj)/c is the propagation time between the 
nodes and /in, 1 (t ) = fi n (t) is a vacuum white-noise 
operator. To simplify notation, we assume quantities 
without index to be the same for all nodes and define 
rotated and retarded photonic operators, i.e., Cj(t) — > 



Ci(t - Zi/c) and similarly for / ini ;, / ou t,i, £%■ All 



a iui L Zi/c 

other quantities are redefined as, e.g., bi(t) — > bi(t — Zi/c) 
and we can then shift the time in the QLEs for node i ac- 
cording to t — > t + Zi/c. Including local dissipative effects 
as above, the full OM QLEs for A — > read 



h{t) = i[Hl ode ,bi\ - 
Cj(t) = i[Hl ode ,a\ - 



L bi - x/7m6(*) . 
/ 2K//in,i(i) 





where the cascaded coupling is encoded in the input- 
output relation 



fin,i+l(t) — fout,i(t) 



(28) 



Now, we can treat the classical laser drives by making 
the replacements c, — > c, + a; and bi — !• bi + fy, and 
from the requirement of vanishing classical forces on the 
shifted operators we obtain the system 



= - (i A c + n)a t - ig (ft + ft* ) a t + £f , (29a) 



ft = -(iu r + 7 m /2)ft - ig \otiY 



(29b) 



2k f y^yj^j ot-i is the effective coherent 



Here, £f = £ 

drive for node i, which is the superposition of the lo- 
cal drive and the output of previous cavities. By taking 
this accumulation into account, it is in principle pos- 
sible to control the classical fields in the cavities inde- 
pendently. However, by considering schemes with more 
than one cavity per node (42J [43], one could also pre- 
vent the classical control signals from propagating down 
the fiber in the first place and thereby achieve indepen- 
dent tunability more easily. The steady state for each on 
is given by Eq. (11) with the replacement £ — > £f s and 



similar statements on adiabaticity apply. Provided that 
\ai\ 2 <§; 1, we may linearize the OM coupling as before 
and the resonators and cavities are then described by the 
linear system 



in Eq. (151. The vector I l (t) 



contains the fiber-input of node i, while R 4 (t) = 

(■/ym&> V / 2«o/o,i, ^l^S,}, VZ^ofliV summarizes local 
noise inputs. 



A. Ideal effective qubit network 

We derive an effective description of the qubits on 
the basis of the assumption that their coupling to the 
resonators is slow compared to the OM dynamics (cf. 
Sec.[H]). This allows us to eliminate the OM degrees of 
freedom and most importantly, we expect the qubits 
to inherit the cascaded nature of the coupling present 
between the cavities. We have seen before that a sin- 
gle OMT mediates a qubit-fiber coupling ~ yT, and 
we hence expect the effective cascaded coupling between 
qubits i and j due to emission and re-absorption of a pho- 
ton to scale as y/TiTj. The general elimination procedure 
is presented in App.[S] and [B] and for later convenience 
in deriving the state transfer protocol and discussing im- 
perfections, we present the results in a master equation 
(ME) formulation for the reduced qubit density operator 

To focus on the key points, we postpone the general 
situation to the next subsection and first discuss the ide- 
alized case in which we (i) take the OM coupling in RWA 
(C = 0) and (ii) assume all additional decay channels and 
noise sources to be zero. The general effective ME (BI) 
then assumes a very simple form: 



fi — £qM + AdealA* • 



(31) 



Here, £ M = — i Si ^q[ a h lAI^ 1S * ne renormalized free 
qubit Liouvillian with uj q — u> q + Ao,i, where the shifts 
A i arc defined below. The second term describes the 



ideal cascaded interaction 



and reads 



H c ff = 



(32) 



(33) 



where S = S; y/Tia~ is the collective jump operator and 
the single-qubit decay rates I\ are given below. The first, 
anti-hermitian term in the effective Hamiltonian H e s en- 
sures that Adeai is of Lindblad form with a single jump 
operator S and the second, hermitian term describes the 
coherent part of the fiber-induced dynamics. In writing 
the above equations, we have further absorbed phases 
e~' l9i a~[~ —> <j~ into the qubit operators to simplify nota- 
tion |^] Note that the two-qubit terms in the hermitian 



v l (t) = - MV(t) - K l (t) - V2^V(t) ■ 



(30) 



The two-qubit terms are determined by the cascaded coupling 
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and non-hermitian parts of H c ff have the same magnitude 
and a specific phase-relation, such that they actually in- 
terfere to produce the cascaded coupling. This is more 
evident when rewriting H c h as 



H. 



eff 



E 



r, 



3 °i a j 



(34) 



*>3 



Apart from generating single-qubit decays, this Hamil- 
tonian may transfer excitations from qubit j to a sub- 
sequent qubit i > j, but not in the other direction. As 
expected, these cascaded interactions take place on the 
scale y/TiTj and we will exploit them to perform a state 
transfer below. 

The effective quantities entering the above ME are 
given by the dynamical properties of the underlying 
multinode OM system. To determine the cascaded dy- 
namics we use the general relations from App.[B"| 

A 2 A 2 

I\ + 2iA , s : = —X u {u}q) , Jij = —Xij((jj q ) , (35) 

where the cascaded coupling J i3 determines the two-qubit 
terms in H c ff. The correlation function Xij{ui) is gener- 
ally defined as 

/>oo 

X ij (u)= dr([6 l (r),6V0)]) free e iwT ) (36) 



where the subscript "free" denotes the expectation value 
in the steady state of the OM dynamics in the absence 
of the qubits as described by Eqs. (28) and (30). It can 



be evaluated using the quantum regression theorem [49] 
with the initial condition ([&i(0), &J-(0)])f rcc = 5y. Defin- 
ing the OM response matrix of node i as A l (uj q ) — 
(M* — iujt)^ 1 , we find in general that Xu(u) = A^u) 
and also that Jij — for i < j, which yields the uni- 
directionality of the coupling. In contrast, for i > j the 
correlation function describes the propagation of a signal 
from node j to node i including various filtering and ab- 
sorption effects taking place at intermediate nodes (see 
App.[B] for details). For the idealized case considered in 
this section, these effects only amount to a phase and we 
obtain |Jy| = ^JTiTj as expected. 



B. Full master equation 

The ideal picture presented in the previous subsection 
has to be refined in order to discuss the impact of the 
various imperfections. As has already been discussed 



Jij introduced in Eq. j35[ l. In OM RWA and for -y m = ko = it 
can be written as Jij = yTiTj cxp[i(#; — 0j)] for i > j, where 

Si = <f>i + 5Zn=i 2<£n wr th 4> n = a<rg{'^2l(' J 3)} ant ^ ^« taken to 
be real for simplicity. This structure allows to absorb e~ lBi a^ — > 



If 



(=) 



FIG. 3. (Color online) Accumulation of noise in a cas- 
caded network where. Thick, wavy arrows denote noise pho- 
tons (red, single-headed) and noise phonons (green, double- 
headed). Each node emits noise photons generated by 
non-RWA scattering events or by up-conversion of thermal 
phonons. At successive nodes, these noise photons may be 
down-converted and lead to decoherence of the qubits (thin 
wavy arrows). 



in Sec|ITJ the effects of the thermal noise and non-RWA 
corrections for a single node are two-fold: They lead to 
heating of the attached qubit as well as contamination of 
the OMT's output with noise photons. In a multiqubit 
network, these noise photons naturally affect successive 
nodes and when moving down the fiber we expect an ac- 
cumulation of noise, as illustrated in Fig. [3] The general 
ME for an iV-node network capturing all these imperfec- 
tions is derived in App.[B]and reads: 

A = \ E { - ^H.m] + TiV[ar}^ (37) 

i 

i>j 

+ lj2TiNi([[at,v],*-] + [at,[»,a7]]) 



E 



Di 



> Ml 



where 2?[a]/i = 2afia^ — a' a fi — /i a^a is a Lindblad term 
with jump operator a. The first two lines contain the 
fiber-induced decay and cascaded coupling of the qubits 
as discussed in the previous subsection, with the associ- 



ated rates Ti and Jij given in Eq. (35). The renormal- 



ized qubit frequencies are given by uj q = uj q + Aj, with 
A, = Ao,i + Ath,jj where the second, thermal contribu- 
tion A t h i is negligible for our purposes (see the discussion 
in App[B[). In contrast, the last two lines describe addi- 
tional decoherence processes, characterized by effective 
bath occupation numbers Ni and the rates of correlated 
diffusion D™. These quantities are given by 



TiNi = ^Y lt {uj q ), 



D; 



A 2 



YM, (38) 



where (w) is determined by the OM dynamics accord- 
ing to 



Yijiu) = 



■<6l(r)6 i (0)) fa 



-itUT 



(39) 
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and the property Yij(uS) — Y^ui) ensures that the Ni are 
real. Finally, we note that the explicit Lindblad form of 



Eq. ( 37 1 is in general not very illuminating and also te- 



dious to compute except for special cases as given below 
or in the previous subsection. Nevertheless, it is ensured 
by the property {b\{r)b 3 ;(0)) free = (^(0)fe J (-r)) frcc of the 
OM correlation functions that the qubit density opera- 
tor (i remains positive semi-definite under the evolution 
described by Eq. @ (see Ref. [54]). 

In the previous subsection, we assumed the idealized 
conditions of vanishing intrinsic decays kq — j m = 0, 
and C = 0. In this limit, we obtain Ni = 0, = 0, 
and | Ja | = yjYiYj, which allows one to rewrite the result 



Eq. (37 1 in the simple form of Eq. (31 1. We will now relax 



these restrictions to discuss the influence of the various 
imperfections and give the conditions under which their 
effects are small. 



1. Intrinsic cavity decays 

For non- vanishing intrinsic cavity decay kq only a por- 
tion r\ = Kf /k of the qubit excitation is actually emitted 
into the fiber, while the rest gets lost to other channels. 
For the simple case of two nodes and ( = 0, the ME 
can be written exactly as the sum of a reduced cascaded 
interaction plus additional on-site decays: 

ft = £ A* + ??Ad e al/U + (1 - T}) -^^i ]M > ( 40 ) 



with given in Eq. (20). For more than two nodes the 



ME is less simple, since every cavity provides an ad- 
ditional decay channel with an associated jump oper- 
ator involving all previous qubits. Therefore, we only 
note that the cascaded couplings decay exponentially 
with the number of intermediate nodes, i.e., \Jij\ = 

yTi ■ ■ - tj+i\ x/Fj where ti(aj q ) = C^i^q) are trans- 
fer amplitudes of the intermediate cavities with the ma- 
trices C 1 {uj) defined in App.|B| For kq <C Kf and A c = u) r 
we find 1 — < 2k /k/. 



2. Non-RWA corrections 

The most important effect of the non-RWA terms in 
the OM coupling is the appearance of additional con- 
tributions to the single- and multiqubit noise terms due 
to Stokes scattering, as has already been discussed for 
the single-qubit case in Sec.|ll] (see Eq. (22)). Since ex- 
plicit expressions for these contributions are lengthy and 
not very illuminating, we only mention that their effect 
scales as k 2 /uj^, |G| 2 /w^ and they are therefore small in 
the resolved-sideband limit. Note, however, that the non- 
RWA terms in the various nodes are in principle phase- 
coherent, which can lead to non-trivial interference ef- 
fects. Apart from the additional noise, one also obtains 



many rather quantitative corrections of the decay rates 
and cascaded couplings as compared to the RWA case, 
such as the fact that the relation | | ps yTVEj becomes 
approximate. However, the general picture of emission 
and reabsorption of qubit excitations as conveyed by the 
ideal cascaded Liouvillian Adcai remains valid. 



Accumulation of noise 



In Sec. |II E| we have discussed the photonic noise 
present in the output of a single OMT. Clearly, in the 
current multinode setting these noise photons will prop- 
agate down the fiber and affect successive qubits (see 
Fig. [3]) . Since the local noise sources (encapsulated in the 
Rj(t) in Eq. (30)) are uncorrelated, we expect the effec- 
tive bath occupation numbers to scale with the number 
of preceding nodes, i.e., 2Vj cx i. To be more specific, we 
approximate them as 



Ni « N Q , + N c 



(41) 



Here, No,i is the local contribution as given before in 



Eq. ( 22 ) and the effect of previous nodes is summarized 
by the cascaded occupation number N c i . Neglecting a 
non-RWA pre- factor to N Ci i, it is defined as the spectrum 
of the cavity input evaluated at the qubit frequency, i.e., 
N c ,i = JZo dT (fLA T )f^(°))irc C e-^ T and we have 
iV Cj i = for the first node. 

More explicitly, for the case of N — 2 nodes, the cas- 
caded noise is simply determined by the output of the 
first cavity and it can be expressed as N Ci 2 = 2kj C\\{u} q ), 
with the steady-state correlation matrix C 11 (w) defined 
in App.|5] Assuming that the first qubit is on resonance 
with one of its associated OM normal modes, we obtain 
for k, \G\ < uj r = A c : 



11 

K 



, 7m N„ 

Top 



7op 



(42) 



where x = 1 f° r weak OM coupling (\G\ < k/2) and 



i|G|" 



l(\G\ 



16' 



") < 2 for strong OM coupling 



(\G\ > k/2). When comparing this to the naive estimate 
for the excess noise given in Eq. ( 24 ) , we note that we 
have obtained the scaling N Cj 2 



0,1- 



This means that the OMT at the second node provides 
exactly the filtering needed to suppress the noise origi- 
nating from the first node (cf. Sec. HE). Therefore, we 



conclude that exchanging a single photon between two 
nodes of the quantum network is in principle possible. 

For completeness, we briefly consider the general case 
of many nodes. To get a quantitative picture, it is again 
instructive to consider the simple case of £ = and 
Ko = 0, such that for j m — > the only imperfection is 
mechanical diffusion at rate jmNm. In this case, Eq. (41 ) 
holds exactly and we obtain 



=X>/i^?i(" 



7, 



A" 



(43) 
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FIG. 4. (Color online) State transfer in a two-node cas- 
caded quantum network, (a) Schematic illustration of time- 
symmetric photon wave-packet (see text). (b),(c) Exem- 
plary pulses Yi(t) given in Eq. ( |50[ ) and Eq. ( |51[ l, respec- 
tively, and resulting occupations puTj 2 obtained numerically. 
In both cases, the parameters have been adjusted such that 

Kfe)! 2 < io- 2 . 



where the appearance of the response matrix element 
■^21 reflects the fact that the thermal noise is first 
filtered by the OMT at node n before it propagates 
down the fiber to affect node i. On resonance, the 
terms in the sum scale as im^m/lop an d we thus ob- 
tain the expected accumulation of thermal noise. For 
identical nodes, this means Nj ~ i "/ m N m and a similar 
statement holds for the multiqubit diffusion terms, i.e., 
Dij r~j min(i, j) 7 m -/V m r/7 op . We finally note that within 



the present approximations, all of these results can be 
expressed compactly by rewriting the last two lines of 
Eq. (37 1 as a sum over the thermal noise sources: 



N n ( \ l"S t 

i iy m V L L 



PI, [M, S n]]) 



This expression is explicitly in Lindblad form and the 
jump operator associated with the thermal noise origi- 
nating from node n involves all subsequent qubits accord- 
ing to S„ = | J2i>nO~ii( UJ q))* a i ■ Here, the matrices 
T m (u/) describe propagation from node n to node i (see 
App.rBI) and we obtain ~ if the qubits are 

resonant with the OM normal modes. 



C. State transfer protocol 

Up to now, the dynamics of the cascaded qubit network 
has been discussed on general grounds. We now turn to 
a specific application, namely the transfer of a quantum 
state from one qubit to the next. This has been the 
topic of our recent work |42j and we provide here the 
details of the protocol. The main problem in transferring 
a quantum state |i/>o) = ot\Q) + /3|1) between two nodes 
of a cascaded network according to 



|^o>i|0) 2 -> |0)i|Vo> 



(44) 



lies in the possibility that the photon emitted by the 
first node could pass the second node instead of being 
reabsorbed (see FigQa)). It was first realized by Cirac 
et aL[Tl] in the context of atomic cavity QED that such 
photon loss can be avoided by choosing appropriate time- 
dependent control pulses in the nodes, which leads to 
a deterministic state transfer protocol. In our setting, 
which is closely related, we can adiabatically tune the 
effective qubit-fiber couplings Ti(t) as well as the bare 
qubit frequencies We proceed by first deriving the 

ideal pulse-shapes needed for the state transfer and then 
comment on their realization. 



1. Pulse Shapes for state transfer 

We ignore all imperfections for the moment and base 
our derivation on the ideal cascaded ME as given in 
Eq. ^ for N = 2 qubits. Along the lines of Ref. P] 
we require that the system remains in a pure state 
fi(t) = \ip(t))(ip(t)\ during the entire evolution, which is 
equivalent to requiring the output of the photo-detector 
indicated in Fi g.[4^ a) to be exactly zero. It is clear 
from Eqs. < 32 ) , (|33[) that this is the case if we can en- 



force the so-called dark state condition S(t)\ip(t)) = 
for all times, where the time-dependence of the jump 
operator is attributed to the time-dependence of Ti(t). 
Under this constraint the evolution of the system is 
completely characterized by the Schroedinger equation 
d t \i)) = -i{H eS + H Q )\i)), with H eS given in Eq. ([33]) 



and Hq = J2i & l q {t) a l/2 containing the (renormalized) 
qubit level splittings. We further re-introduce the phases 
9i that had been absorbed in Sec. |III A\ and expand the 
state vector in terms of three time-dependent amplitudes 
u,v%,V2 as follows: 

|V(*)> =cm(i)e l *+ (t) |00) (45) 
+ p L (t) e^- {t) |10) + «a(t) e »-(*)+i*(t) |oi)l . (46) 

Here, the phase factors have been introduced for later 
convenience according to $± (t) — | ^ ds (uj^(t)±£jq(t)) 
and 4>(t) = arg{J 2 i} — 62 — 6\. The dark-state condi- 
tion S\%jj) = now reads vTT w i + vT^t^ = and the 
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Schrodinger equation amounts to u = and 



-vi , v 2 



£2 
2 



iS(t)j v 2 - v /r i r 2«i , 



(47) 



where 5(t) = t~j q (t) — u q (t) + 4>(t) is an effective detuning. 
We will first derive the pulses for S(t) = 0, which is a 
generalized resonance condition taking the varying phase 
of the cascaded coupling into account. 

Achieving perfect state transfer means that we find 
time-dependent r, (t) such that the solution of Eq. ( 47 ) 



satisfies the dark state condition as well as the boundary 
conditions 

Vi{U) = \v 2 {t f )\ = l, v 1 (t f )=v 2 (U) = 0, (48) 

where ti and tf are the initial and final times, respec- 
tively. In reality, however, one will have to tolerate slight 
violations of these boundary conditions due to finite pulse 
lengths, etc. To find suitable pulses we make use of the 
following time-symmetry argument |lll 155] : If a pho- 
ton is emitted by the first qubit, then, upon reversing 
the direction of time, we would see a perfect reabsorp- 
tion. We can exploit this by ensuring that the emitted 
photon shape is invariant under time-reversal and use a 
time-reversed control pulse for the second qubit. As a 
result, the absorption process in the second node is a 
time-reversed copy of the emission in the first and may 
thus be - at least in principle - perfect. In the present 
formulation the role of the photon shape is played by 
a(t) = \/T\vi and we can generate a class of pulses by 
requiring its time-derivative to be of the form d = f(t)a 
with some f(t) satisfying f(t) = —f(—t), provided we 
assume tf = — i,. Upon using Eq. (47) this yields a dif- 



ferential equation for the pulse shape: 



r a = r 2 + 2f(t)r 1 



(49) 



After choosing a suitable (i.e. bounded and positive) so- 
lution one may use it to calculate v\ (t) and the quantities 
in the second node are then given by v 2 (t) — —Vi(—t) 
and r 2 (i) = Tx(—t), which solve the Schrodinger equa- 
tion and satisfy the dark-state-condition. Note also, that 
related time-reversal arguments can be employed for re- 
absorbing almost arbitrary photon wave-packets, as has 
been discussed for formally similar atomic systems |56j . 

We assume that there is a maximal achievable decay 
rate r max and give two example pulse-shapes. First, for 
f(t) = — T max sign(t)/2 the wave-packet has the shape 
a(t) cx e~ rmax l*l/ 2 and a solution of Eq. (49) is given by 



r 1 (t<o) = r ma 6 r 1 (<>o) = r„ 

2 — e 1 maxI 



r ma xt 



(50) 



The resulting dynamics is illustrated in Fig.[4](b), where 
the pulse-length has been chosen to be T p = tf — ti = 



8/r max such that |fi(i/)| 2 < 10 2 . As a second example, 
we consider f{t) — —ct for c > and obtain the solution 



ri(t) = 



exp(— ct 2 ) 



Ti(0) 



(51) 



where the parameters c and Ti(0) have to be adjusted 
such that the boundary conditions (48) are met. Again, 
we display an example with |i>i(t/)| < 10~ 2 in Fig.Ulc). 
Note that the pulse-length T p — 12/r max is slightly 
longer than for the previous pulse, since the maximal 
decay is reached only shortly. 



2. State transfer between two nodes 

To discuss the implementation of the state transfer pro- 
tocol in more detail, we focus on the pulse shape given 
in Eq. (51). We expect the other pulse (cf. Eq. (50)) to 



perform worse in terms of non-RWA noise due to the fact 
that the maximum decay is required over a long period. 
The following discussion is, however, valid in general. 

The task is to tune the OM parameters in such a way 
that the resulting effective decays Ti(t) in the two nodes 
follow the prescribed time-evolution Ti(t), and we first 
consider varying the OM couplings Gi(t) as depicted in 
Fig.[2jc). The key step is then to solve the equation 
Ti(Gi(t)) — Ti(t) with k = 7, n = for the desired 
temporal variation of the OM coupling Gi (t) , where we 
use the full relation Eq. ( 19 ) to also capture non-RWA 
corrections. Note that the renormalization of A c has to 
be taken into account and we have further ignored the ef- 
fective qubit shifts as well as the varying phase of the 
cascaded coupling. This means that the resonance con- 
dition is in general not satisfied, i.e., we have 8(t) 7^ 0. 
In order to correct this, we can adjust the bare qubit fre- 
quencies and in a first iteration set Lo l q — > uj' l q + 6uj' q (t), 
with 5u\{t) = Suj 2 (t) = S(t)/2. This shift is of order 
•^ 2 /7op and the new uj q approximately satisfy the reso- 
nance condition with corrections being of higher order in 
A. In principle, shifting the bare qubit frequencies also 
modifies the Ti(t) which therefore deviate from the pre- 
scribed Ti(t). However, this deviation is negligible since 
the shifts dui q are small compared to the width 7 op of 
the resonances in T(u> q , G) [^] We have thus obtained the 
pulse (Su} q (t), Gi(i), Sco 2 (t) : G 2 (t)) that realizes the de- 
sired control pulse I\(i) and satisfies the resonance con- 
dition. 

An example for a resulting pulse is displayed in 
Fig.[5ja) and we note that the adiabaticity conditions 
are well satisfied: Since G(t) is tuned from up to its 



3 To satisfy all of these constraints exactly, one may perform a 
minimization of the functional C[<5u>*(t), Gj(t)] = f t f dt [(<5(t)) 2 + 
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maximum value in about one quarter of the pulse-length 
T p , we have on average G/Gm&x r ^ > X/Tp oc r max 

< 7o P , 

as is required for the adiabatic elimination to hold and 
is therefore valid by assumption. In addition, for trans- 
lating the desired pulse G(t) into applied laser power via 
Eq. (11 1, we need G/G m & x <C k, which is implied by the 
previous relation. To assess the impact of imperfections 
we can now evaluate all quantities in the effective ME 
(371 for potentially finite j m , kq, N m , as will be done 



below. Finally, as an alternative to tuning the OM cou- 
plings Gi , we can also vary the cavity frequencies in order 
to regulate the effective decays, as depicted in Fig. 2^d) 
and a corresponding pulse A C) j(t) is shown in Fig. [5 b). 



3. Discussion 

The quality of the state transfer can be estimated by 
averaging the fidelity 



(52) 



over all transmitted single-qubit states k/>o) = a|0)+/3|l). 
Here, n(i) is the solution of the ME (37) with initial con- 



dition /i(0) = \ipo)(ipo\ ® |0> <0| and \ip ) = ae 4 *+ (t) |0) 
^.f- (*)+««*) |1) is the target state including the various 
phases that are accumulated due to the varying qubit 
frequencies and phase of the cascaded coupling. We also 
take into account a decay of the qubit coherences ac- 
cording to e~'/ T2 by adding a term ^ V[al]ii to the 
ME, where T2 is the qubit's intrinsic coherence time. The 
leading contributions to the fidelity are expected to scale 



T k, 1 



2 Kg 

3 K 



C 2 



c. 



A 2 T 9 



(53) 



where the Ci are numerical constants of order one and 
the four terms correspond to intrinsic cavity loss, ther- 
mal noise, Stokes scattering (we chose G max ~ k for the 
pulses), and intrinsic qubit dephasing. The constants Ci 
depend on the pulse-shape and can be optimized depend- 
ing on the experimental parameters. Figs.[5]jc),(e) dis- 
play the fidelities obtained in Ref. [32] for varying Gi(t), 
while Figs.[5]jd),(f) show results for varying A Cj j(t). It 
can be seen that the two control schemes do not differ 
significantly, although we will use the former in our esti- 
mates below. 

To conclude this section, we briefly discuss the results 
for an implementation with spin- and charge qubits, the 
details of which are described Sec.|V| For spin qubits, 
we have X/2tt on the order of 50 kHz if we choose 
uj r /2'K « 5 MHz and we assume n/2ir m 1 MHz to be 
in the resolved side-band regime. With expected deco- 
herence and decay rates (Ko,rm)/27r = (50,10) kHz, as 
well as T2 ~ 10 ms, we obtain state transfer fidelities 
of T rs 0.85, where all imperfections contribute approx- 
imately equally. However, with recent improvements of 
T2-times for spin qubits [57], fidelities beyond T ss 0.9 




0.04 K /A 2 T 2 



FIG. 5. (Color online) Performance of state transfer proto- 
col using the pulse of Fig.[4|c). The left and right columns 
show realizations via tuning of OM couplings and cavity fre- 
quencies, respectively. (a),(b) Resulting control pulses for 
k = 0.05uv (see also Fig.[2}. (c),(d) Fidelities as functions 
of k and r m , where ko = I/T2 = 0. (e),(f) Fidelities as func- 
tions of k and I/T2, with kq = F m = 0. For the left column 
we deduce the parameters of Eq. (53 1 to be Ci ~ 4, C2 ~ 1.4 
and C3 ~ 7.5, and for the right column C2 remains the same 
while the others are slightly worse, i.e., Ci « 5 and C3 « 8. 



seem feasible for our parameters. In contrast, the fidelity 
for charge qubits is clearly dominated by intrinsic qubit 
dephasing. In this case, one can choose a larger resonator 
frequency cj r /27r w 50 MHz and k/2it - \/2ir - 5 MHz, 
whereas present-day technology allows for qubit dephas- 
ing times around T 2 ~ 2 fis. Again, we obtain T sa 0.85, 
but here, dephasing alone already accounts for an in- 
fidelity of ~ 0.1. In this case, strategies to overcome 
this problem would have a significant impact on the per- 
formance of the state transfer. In particular, one could 
numerically optimize the pulse-shapes for gate-speed or 
depart from the adiabatic elimination and treat the full 
dynamics of the OMT, yielding faster gate-sequences. We 
conclude that state transfer operations with both types 
of solid-state qubits are feasible with present-day tech- 
nology. 



IV. OPTICAL ON-CHIP NETWORKS 

So far we have considered networks where quantum in- 
formation is transmitted through a continuum of modes 
supported by an optical fiber using pulse-shaping tech- 
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niques. While this is a natural setting for long-distance 
quantum communication, one may also imagine situa- 
tions where the different nodes of the quantum network 
are closely spaced and could instead be linked by optical 
resonators which exhibit a discrete set of modes. This 
scenario could be relevant for quantum communication 
and entanglement distribution within quantum process- 
ing architectures built on a chip, with optical channels 
providing a fast and robust alternative to other communi- 
cation schemes, as discussed in the introduction. There- 
fore, we consider a setup as shown in Fig.[6]Ja), where 
each qubit is linked via an OMT to a "bus" -cavity of 
length L, which may, e.g., be implemented by a short 
optical fiber. Similar to the long-distance network con- 
sidered above, the system is described by a Hamiltonian 
of the form 



H 



tdc 



m 



-link) 



H 



link " 



(54) 



where the structure of H^ odc is the same as in Eq. ([2]), 
with the exception that we do not apply laser drives to 
the nodes (see below). The bus is modeled by an ad- 
ditional cavity with discrete modes separated by a free 
spectral range of Slu/2tt > 1 GHz, which corresponds to 
L < 10 cm. We take this spacing to be much larger than 
the other couplings and decays and therefore, we can re- 
strict our discussion to a single fiber mode with destruc- 
tion operator do and frequency ujq close to the frequencies 
bf c of the nodal cavities. In this case, we obtain 

tfiink - ^o4 d o + i(Se- iULt 4 - £*e iULt d ) , (55) 



h: 



cav-link 



= hi d\ci + h* doc] , 



(56) 



where we have also included a laser of frequency ujl which 
coherently drives the fiber mode. Here, we have neglected 
additional degenerate modes in the nodes which do not 
couple to the fiber. They will not be classically occupied 
and will thus not couple significantly to the mechanical 
resonance (see Sec. II C). 

To proceed, we can follow the basic steps described in 



Sec. Ill First, we introduce the QLEs for the OM degrees 
of freedom to describe dissipative effects [49] : 



d 

Ci 

b, 



-i[do,H'} - KQfd - a/2/«o//o(£) 
-i[ci, H'] - K Ci - \/^Ko/o,i(*) i 




(57a) 
(57b) 
(57c) 
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FIG. 6. (Color online) Phonon-photon bus realized via a fiber 
cavity, (a) Schematic illustration of the setup. Each circle 
represents an OM mode and the long box indicates the fiber 
cavity which is driven by a laser of frequency ujl- Wavy arrows 
indicate cavity decay channels (red, single-headed) and me- 
chanical diffusion (green, double-headed), (b) Level scheme 
for N — 2 identical nodes (black levels) corresponding to the 
results presented in Tab.[l] The two outer levels signify the 
two qubits to be coupled. For N > 2 the modes labeled (a, ±) 
become N — 1-fold degenerate (indicated in gray). The nor- 
mal modes in RWA and for vanishing dissipation are given 



by w± = (pa ± c a )/\f2, as well as Wq = ^-b s - 
u4 = -^(f b s ±~c s + f co), where 5 = VG 2 + K*. 



-co and 



QLEs to vanish then yields (i > 0) 



a 



-(iA cQ + K f)a Q -i^^hioii + £ , 



(58a) 



Here, H' = H — H cnv and k , j m /2 are the decay rates 
of the cavities and resonators with corresponding noise 
operators fo,i(t), 6(0 as introduced after Eqs. ([7]),(|9|. 
In addition, we have introduced a decay rate Ko/ and 
vacuum noise operator fo(t) for the fiber mode. Sec- 
ond, we determine the classical steady state by trans- 
forming all photonic operators to a frame rotating at 
the laser frequency and subsequently making the replace- 
ments d a — s- do + ceo, ^ — >• Cj + cti and 6j — > hi + ft. De- 
manding the c-number contributions to the transformed 



(iAci + no)ai - ih*a - igo&i (ft + ft*) , (58b) 

(58c) 



Q-i = 

ft = -(iw r + 7 m /2)ft - ig \a { \ 



where we have introduced the detunings of the cavities 
from the laser drive according to A c o = ojo — wj, and 
A C i — u>l — uj L . Below, we are interested in the case 
of identical nodes, where the above system admits for a 
uniform steady-state solution with identical a.j and ft. 
For \an\ J > 1 we may now carry out the linearization of 
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the QLEs (57), which amounts to the replacement 



H' — > H n 



E 



2 



crib, 



(59) 



where H osc is given by 



H osc = A c0 dld + Y (4 ^ c * + h - c -J ( 60 ) 

i 

+ [u r b\bi + A cl c\c % + (G,c| + GJci)(6j + b\) . 

i 

Here, G{ — goai is the enhanced OM coupling and we 
have redefined A C i — 2\Gi\ 2 /uj r — > A c j. In summary, we 
have thus obtained a system of qubits which is connected 
by a finite linear network of harmonic oscillators. 

The linear oscillator network H osc can in principle be 
used to resonantly propagate excitations from one qubit 
to another. However, this would require a certain degree 
of control over the various couplings to route an excita- 
tion to its target node and also has the drawback that 
the oscillators' decay channels have a strong effect. For 
these reasons we will focus on an off-resonant scheme: 
The qubits chosen to take part in the gate operation are 
tuned close to one of the normal modes of the intermedi- 
ate OM system described by -ff sc- If A is the detuning 
from that mode they experience an effective interaction 
at rate J ~ A 2 / A, provided A <C |A|. On the other 
hand, a decay 7 of the chosen normal mode will lead to 
an induced qubit-decay Y ~ A 2 7/A 2 , which can be sup- 
pressed for larger detunings according to Y/J ~ 7/A, as 
long as the intrinsic noise-processes of the qubit do not 
play a role (see below). To exploit this scaling, we need 
of course well-resolved OM normal modes. 

A. Two node network 

We illustrate the resulting interactions for N = 2 
nodes and postpone matters regarding the scalability to 
Sec. |IVD) To begin with, we analyze the mode-structure 
of the OM system described by iT sc, which will provide 
the basis for understanding the induced qubit dynam- 
ics. Assuming identical nodes, we introduce symmet- 
ric and anti-symmetric combinations of all system and 
noise operators (i.e. c s / a = (ci ± c 2 )/V2, etc.). The 



second term in Eq. (|60|) then simply reads Kdlc s + h.c 



with K = \/2h and the QLEs for the oscillator network 
(Eq. ( [571 with H' — > H osc ) decouple into a symmetric 
and an anti-symmetric set: 



v"(t) = -M"v w (f) - R"(t) v = 



(61) 



Here, the degrees of freedom have been grouped into 
v s = (b s ,c s ,d ,bl,cl,dl) T and v a = (b a , c a , b\, c\) T , and 
the noise vectors R"(i) are mutually uncorrelated. Ex- 
plicit expressions for the M" and R"(i) can be found in 
App.lO and we take G, K > from now on since their 



phases can be absorbed into the cavity operators. The 
anti-symmetric set behaves like a normal OM system 
as discussed in Sec. [II] and for the relevant strong cou- 
pling case G 3> Ko,7 m , mechanical and optical modes 
hybridize to produce normal modes at u)± ~ uj r ± G. 
Here and in what follows we assume A c = A c o = cu r and 
the results in the present subsection are obtained for OM 
RWA (£ = 0). In contrast, the symmetric subsystem is 
an augmented OM system including also the fiber cavity 
with normal modes at loq = u r and oj± ~ ui r ± 5, where 
5 2 = G 2 + K 2 . Therefore, symmetric and anti-symmetric 
normal modes are staggered as shown in Fig.[6]jb) and for 
G,K ^> Ko,Ko/j7m the modes are well resolved, as can 
be seen from the expressions for their widths 7 listed in 
Tab.U 

We proceed to eliminate the OM degrees of freedom on 
the basis of the assumption that they are detuned from 
the qubits by much more than A. In terms of the new OM 
modes, the qubit-resonator interaction in Eq. (59) reads 



i^int — 



E 

y—s.a 



(62) 



and since the b v are mutually uncorrelated they may be 
eliminated independently, as is done in App.[C] However, 
for well-resolved normal modes it is reasonable to further 
assume that also correlations between modes of the same 
set are small. To be specific, we introduce normal modes 
w 1 ' via v" = U v ~w v ', such that the transformations U v 
diagonalize M v . This yields 



E( 

k 



(Unj k L K , (63) 



where ico" + 7J are the eigenvalues of M u . In addition, 
the qubit-resonator interaction becomes 



' '\» + h.c. 



\» = \{V v ) Xi , (64) 



where we have introduced the effective couplings AJ and 
the primed sum is restricted to normal modes with fre- 
quencies ojj > 0. In the present RWA this is exact, but it 
presents an approximation otherwise. Now, correlations 
between the normal modes of a given set v arise due to 
the noise terms in Eq. (63), which scale as y/"f m , •v/2acq/ ) 



\/2kq. Since these quantities have to be small anyway, we 
neglect these correlations and the normal modes Wj then 
evolve independently. As a result, the general elimination 
procedure as outlined in App.[X]becomes very simple and 
calculating the necessary correlation functions ( |A5[ )-(A8) 
for each yields a ME that is just a sum over the five 
normal modes listed in Tab.[TJ 

fi = -i[j((T 1 + a 2 _ +oW + ),v\ (65) 
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mode (v,j) 



Nu 



iVnon-RWA 



G 2 

4{G±U r ) 2 



(a,±) 



u r ±G 



7m «0 



j m N„ 
47± 



(8,0) 



A"' 



G J 



+ "^-KO/ 



A 



-A 2 -JmNn 

s 2 M 



^K 2 + ( K 



(8,±) 



±5 



G K A 

4^ 7m + T + 2^ Ko/ 



G 
V2S' 



G ~j m N„ 



4 

9 k 



2S 2 275 



, 01 * 

4(«0/ + KO + g 2 Ko) LO 2 \LU 



TABLE I. Characterization of OM normal modes (8 = VG 2 + A 2 ): The different columns display eigenfrequencies u and 
decay rates 7 as obtained by diagonalizing M", as well as quantities entering the effective ME (651 in the independent-mode 
approximation, i.e., effective couplings A, thermal occupation numbers N t h and non-RWA occupation numbers Wnon-RWA- 
The first four columns were obtained in OM RWA and we have generally made use of Ko,Ko/,7m <C G, A and assumed 
A c = A c o = uj r to simplify the expressions. The last column displays results to first order in the non-RWA terms (see text) 
and we show only leading terms in G 2 /a> 2 , K 2 /uj 2 for j m — > and g = G/K < 1. 



Here, J = J2 V j Jj ^ s tne desired effective coupling and we 
have dropped terms renormalizing the qubit frequencies 
for simplicity, since they are the same for both qubits. 
The couplings and decay rates are directly related to the 
normal modes as described in Tab.[l] and owing to our 
independent-mode approximation they have the form of 
simple Lorentzians (AJ = uj q — wj): 



J v — 



3 



(A 



v\2 



Aj(-l) 



8 ( 7 n 2 + (a-) 2 



(A?) 



( 7 j0 2 + (Ap^ 



(66a) 
(66b) 



Hence, we have obtained the expected scaling Fj/Jj ~ 
7J /Aj when close to a resonance. The alternating sign 
of the Jj is rooted in the fact that eliminating a mode 
of symmetry v yields a shift of the operator o~^o~ oc 
(_ l) 5 >"i((7+ cr~ + erfer^). Note that when tuning the 
qubits somewhere in between two modes, these have op- 
posite symmetry (cf. Fig.[6]Jb)) and their contributions 
to J add up constructively, since the relevant detunings 
AJ also have opposite signs. 

Finally, the influence of the thermal noise on the 
qubits is characterized by the effective occupation num- 
bers Nj ~ J m N m /tf (see Tab.|l]for details). This means 
that the thermal noise is suppressed by an OM cooling 
effect, provided that we have j m N m < kq, Kof- 



B. Non-RWA effects 

The effect of the non-RWA terms is to admix a counter- 
rotating portion to each OM normal mode, giving rise to 
new noise terms as well as quantitative corrections. In 
the case of the full OM coupling the eigenmode trans- 
formation U u is not block diagonal anymore and the 
cigenmodes in Eq. ( 63 ) are driven by counter- rotating 
noise operators. To estimate their effect we can intro- 



duce new noise operators with effective non-zero occupa- 
tion numbers by means of a Bogoliubov transformation 
and then proceed as above and treat the normal modes 



independently. For illustration, the eigenvectors of the 
anti-symmetric system M a are, to lowest order in the 
non-RWA parameter Q and for vanishing dissipation, 



1 

71 



(la 



±C, 



(G 



2V2{G±i 



(Si 



±c. 



(67) 



Calculating, also to first order in (, the matrix (U a )~ 1 
from these eigenvectors then allows to estimate the effec- 
tive occupation number due to Stokes scattering: 



Nl 

± ,non- 



G 2 



RWA 



4(G±w r ) 



(68) 



where we have neglected non-RWA corrections to the 
thermal noise. For the symmetric modes the expressions 
become quite lengthy, so we display in Tab. [I] only the 
leading order terms in G 2 /uj 2 , K 2 /uj 2 . 



C. Discussion: VSWAP gate 



The effective two-qubit interaction obtained in Eq. ( 65 ) 



can be used to swap the states of the qubits or to perform 
the entangling VSWAP operation, which amounts to ex- 
posing the qubits to this interaction for a time t g = 7r/4J. 
To asses the optimal operation point for such operations 
we consider the creation of a maximally entangled state 
\ipo) = (|01) — i|10))/\/2 from the initial product state 
1 1) by applying the VSWAP gate. For this purpose, 
we also account for a finite qubit coherence time T 2 by 



adding the term ^ Si^t !] to the ME (65), which is 
then solved to first order in the Lindblad terms. The re- 
sulting approximation of the fidelity J- — (ipo\fi(t g )\'ipQ} 
for this operation can be written as 



8 ^ \J\ K i ' 8 J T 2 



(69) 



which is a function of the qubit frequency (cf. Eq. (66)). 
Generally speaking, for K > G we expect higher fidelities 
if the qubits are tuned in between the central resonances 
(s,0), (a,±), as can be argued from Tab.[TJ Compared 
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to the outer resonances (s,±), they have larger weights 
A 2 and thus larger J. Further, their width 7 is slightly 
smaller for the relevant case 7 m <C kq rs ko/. 

In the limit where the induced decay T dominates over 
intrinsic dephasing, i.e., T 1/T 2 , we see from the sec- 
ond term in Eq. ( 69 1 that there is a point between each 



pair of neighboring resonances where J/Y is optimal. For 
G = K such points can, e.g., be found half-way be- 
tween resonances (s,0) and (a,±) at u q « uj r ± G/2. 
By taking only the two closest resonances into account 
(which slightly over-estimates the fidelity) we obtain for 

Kq = K f <C G 



To 



1 



7T Ko 

2 G 



1 



7i 



,A„ 



Kq 



2wl 



7T G 

2 AM 



(70) 



In this limit, the infidelity is thus dominated by cavity- 
induced qubit decay (first term in parentheses), provided 
the thermal noise of the resonator is suppressed due to 
lm,N m < kq. The last term in parentheses represents 
non-RWA corrections and is small for moderate cou- 
plings, as required by stability constraints (cf. Sec. [ill). 

In the opposite limit, where intrinsic qubit dephasing 
dominates (1/T 2 3> T), we have to tune very close to a 
single resonance to obtain large J and hence a fast 
gate. There is, however, a trade-off between minimizing 
the dephasing (small A) and the optomechanically in- 
duced decay T (large A). The optimal detuning in this 
case is given by A opt « [(AJ) 2 7 V(1 + 2N»)T 2 /2] 1 / 2 which 
yields a fidelity of 



opt 



/2 7 J(l + 2Ap 
(A>) 2 T 2 



(71) 



There is, however, also the limitation A op t 3> A im- 
posed by the elimination procedure. 

Figure[7]shows the fidelities for an implementation with 
spin and charge qubits as discussed in more detail in 
Sec.[v] Apart from the expression ( 69 1 based on the 



independent-mode approximation (dashed) we also dis- 
play the result of the full elimination, Eq. ( C8 ) , show- 
ing reasonable agreement. For an ideal superconducting 
qubit, the fidelities are rather good and with T 2 w 2^s 
we still obtain T > 0.85 for the parameters of Fig.^a), 
which shows that the performance is limited by intrin- 
sic dephasing. On the other hand, for the spin qubit 
a resonator of lower frequency is better suited (yielding 
higher zero-point motion), which also puts limits on the 
OM coupling G. As a result, the leading term in Eq. ( 70 ) 



proportional to kq/G is larger than for the charge qubit, 
meaning that the intrinsic cavity decays become a more 
severe limitation in this case. For an ideal spin qubit we 
obtain fidelities of ~ 0.9 for parameters as in Fig.[7jb). 
Intrinsic dephasing then has a less severe impact than 
for the charge qubit and for T 2 ps 10 ms we still find 
J- pa 0.85 — 0.9. Strictly speaking, these estimates are 
only valid for small infidelities, but they are neverthe- 
less useful for identifying interesting regimes of operation 
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FIG. 7. (Color online) Fidelities for generating a maximally 
entangled state between N — 2 nodes (see text). Solid curves: 
results of the full elimination presented in App. [C} d ashed 
curves: approximate expressions according to Eq. (691. (a) 
Parameters suitable for charge qubit: uj r /2ir — 50 MHz, 
A/2-7T = 3.5 MHz, G — 0.2a;,., K = 0.25o; r , and decoherence 
rates (k , «o/, r m )/27r = (50, 25, 10) kHz. The blue (dark 
gray) curves correspond to vanishing dephasing I/T2 = 
and the green (light gray) curves to T2 = 2 fjs. (b) Pa- 
rameters for spin qubit: to r /2n — 7.5 MHz, \/2n — 40 kHz, 
G — K = 0.25av and decoherence rates as in (a). Blue (dark 
gray) curves: I/T2 = 0, green (light gray) curves: T2 = 10 ms. 



and we recognize the two limiting cases discussed above. 
As in the case of the long-distance state transfer, one 
might - especially in the case of the charge qubit - profit 
from leaving the adiabatic regime in order to shorten the 
gate-sequence and to reduce the impact of intrinsic qubit 
dephasing. 



D. Outlook: Scalability 

When considering an A^-node network the most 
straight-forward way to perform a gate operation be- 
tween two given qubits is to detune the nodal cavities 
of all other nodes sufficiently far, such that the network 
reduces to the two-node problem, which leads us to ex- 
pect similar fidelities. Since this requires only moderate 
temporal control of the cavity frequencies w*, this could 
be achieved by applying, e.g., strain or heat to those 
cavities not meant to participate (see, e.g., Ref. 58] for 
tuning of the toroidal cavities discussed in Sec.[V). Note 
that if the fiber has to be extended to accommodate the 
new nodes, the cavity-fiber couplings scale as h cx 1/yL, 
which has to be taken into account. 

We now consider A-qubit interactions and the basis 
for the discussion is again the normal-mode structure of 
the underlying A-node OM network. For identical nodes, 
only the "center-of-mass" cavity mode c\ = X)i c */v / A 
couples to the fiber while all other modes decouple. As 
a result, the level scheme in Fig.pfb) remains valid with 
the modification that the modes (a, ±) become (A — 1)- 
fold degenerate (indicated in gray) and we assume that 
the fiber is extended with every node added, meaning 
h cx 1/VN and hence K = \fNh = const. As could 
be expected, the calculation presented in App.[C] shows 
that the induced coherent qubit dynamics is given by 
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FIG. 8. (Color online) Realization with toroidal cavities and 
doubly clamped beams, (a) Illustration of a single toroidal 
cavity coupled to an optical fiber as well as to a doubly 
clamped beam, (b) Coupling of the doubly clamped beam 
to a spin qubit via magnetic fields, (c) Coupling of the dou- 
bly clamped beam to a charge qubit via electrostatic forces. 



JnS + S~ , where S~ = J2i a T an d in general Jn oc 1/N. 
For illustration, we note that if the qubits are tuned close 
to the symmetric mode (s,0) at w r , we can neglect the 
anti-symmetric ones provided G > Nkq and write the 
effective iV-qubit ME (|C2l as 



(i m — i [ JjyS + S , fi\ 

r, 



(72) 



(N coll + l)V[S-}fx- 



^N coll V[S+]fx, 



where the coefficients are related to the two-node re- 
sult via J N ~ 2J/N and r col i ~ 2rg/A, and N coll 
is an effective occupation number. We hence obtain 
Fco\\/Jn ~ const and conclude that decays induced by 
the OMTs do not limit the number of nodes, provided 
we are close enough to the mode (s,0). In contrast, when 
tuning to the anti-symmetric modes, there are always 
N — 1 decay channels present, leading to T / oc N . Fi- 
nally, note that even for the symmetric modes there are 
always the intrinsic decoherence processes of the qubits 
that set an absolute time-scale and limit the number of 
nodes via the requirement J/vT 2 oc JT 2 /N ^> 1. 



V. IMPLEMENTATIONS 

In the preceeding sections, we have mainly discussed 
the OMT on general grounds and conditions for low-noise 
operation have been given. Here, we describe the realiza- 
tion of the OMT within specific physical systems and the 



resulting numbers have already been used in Sec. Ill and 
|IV| to illustrate the performance of the proposed schemes. 
Generally, implementing the OMT poses the challenge of 
combining high-Q optical cavities with high-Q mechan- 
ical resonators exhibiting an appreciable zero-point mo- 
tion, i.e., a low effective mass. Further, the mechanical 
resonators must interact with the qubits, which requires 
suitable coupling schemes. Linking the different nodes 



of a quantum network can be achieved with standard 
optical fibers, while for short distances, integrated nano- 
photonic circuits [131 [SS] might be a promising alterna- 
tive, which we, however, do not explore here. 



A. Toroidal cavities and doubly clamped beams 

As a possible candidate system for the OMT we pro- 
pose to use a setup similar to the one recently demon- 
strated by Anetsberger et al. [60 . There, a doubly 
clamped beam of stressed SiN is positioned in the evanes- 
cent field of a whispering-gallery mode (WGM) sup- 
ported by a microtoroidal cavity made of silica. This 
arrangement gives rise to the standard OM coupling and 
the toroid is further coupled to a tapered fiber for driv- 
ing and interrogation as shown in Fig.[8^a). To obtain 
a resonator-qubit coupling we propose to attach a cou- 
pling element to the resonator (e.g. a magnetic tip or an 
electrode as discussed below), which resides sufficiently 
far outside the evanescent field of the WGM such that 
the optical quality factor is not degraded. To be specific, 
we envision to use a toroid of diameter 2R 20 pm and 
the second harmonic mode of a doubly clamped beam of 
length I w 15 pm positioned as in Fig.[8]ja). The coupling 
to qubit and cavity can then occur at different antin- 
odes and the coupling element is located at a distance of 
d > 2.5 /im from the rim of the toroid, which is sufficient 
for our purposes (see below). 

The elastic modes of a doubly-clamped SiN beam of 
length I, width w and thickness t are well described by 
standard elastic theory [BUI [5T]. For the case that the 
beam is exposed to high tensile stress, the modes re- 
semble those of a vibrating string with frequencies given 
by uj n = &o/ Pm Ttn/l, where n — 1,2,... labels the 
mode, p m is the mass density and <7n the internal ten- 
sile stress. In contrast, for low internal tensile stress and 
w = t the frequencies scale as w„ = c^^E/l2p m t/l 2 
with ci w 4.7300, c 2 ~ 7.8532 [62], where E is Young's 
modulus. We are interested in the second harmonic of 
the beam and hence set uj t — lj2- When quantizing the 
resonator, the effective mass of each mode is given by 
an integral over the mode's displacement profile u n (z) 
according to m e s — ml^ 1 J Q dzu^(z), where m = t 2 lp m 
is the real mass. For the relevant mode we find that 
TO c ff m/2 in both stress regimes. By tuning the in- 
ternal tensile stress, one can interpolate between the two 
regimes and values of p m = 2800 kg/m 2 , E = 160 GPa, 
and <7rj up to ~ 1 GPa have been demonstrated (see 
Refs. (S3 EI])- Therefore, w r /27r « 5 ... 50 MHz for a 
beam with dimensions (l,w,t) — (15,0.05,0.05) /xm is 
realistic and the respective zero-point motions are in the 
range ao = y/H/2m c f[U) r 1.8 .. . 0.6 x 10~ 13 m. Demon- 
strated Q-factors for these devices range from 10 4 to 
1.5 x 10 5 EQlET] and for Q = 2 x 10 5 we obtain, e.g., a 
mechanical decoherence rate of j m N m /2ir s» 10 kHz at a 
support temperature of T ~ 100 mK. 

The microtoroidal cavity is modeled by a pair of 
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(by symmetry) degenerate counter-propagating WGM 
modes, whose field distributions we denote by El^t). 
The theoretical quality factors (including radiation loss 
and absorption in silica) are beyond 10 10 , whereas 
demonstrated values reach up to Q — 4 x 10 s [53"] . 
When the SiN beam is brought into the near-field of 
these modes, they generally experience a frequency shift 
as well as a cross-coupling. Simple perturbation theory 
on the level of the wave equation shows that the rele- 
vant overlap integrals are 1^ oc J* resonator d 3 r AeE*E v 
with Ae = ri 2 cs — 1 and n rcs the refractive index of 
the resonator. However, a resonator oriented as shown 
in Fig.[8](a) is quite different from a point-scatterer and 
hence does not scatter between the modes. To be specific, 
we obtain the scaling Ilr/Irr ~ exp[— k 2 R/x] <C 1, 
where k — 2tt/X c and x _1 = 250 nm is the field decay 
constant outside the cavity [50] , Therefore, we can take 
the OM coupling to be of the type assumed in Eq. ([3| 
with g oc dI RR /dz. 

We now turn to linking the different nodes: When 
the toroidal cavities are side-coupled to the fiber, the 
modes of a given propagation direction experience a uni- 
directional coupling, as needed for the cascaded quantum 
network of Sec. |III| However, for standing wave modes, as 
they occur when the degeneracy of the WGM-modes is 
lifted by an amount h > k, G or alternatively in Fabry- 
Perot or photonic band-gap cavities, the cascaded cou- 
pling has to be achieved by additional non-reciprocal op- 
tical elements, such as optical circulators or Faraday ro- 
tators. While this may complicate the structure of the 
network, the idea of the OMT remains valid. 

In the on-chip setting, the nodes are linked by a cavity 
of length L and refractive index n^, such that its mode- 
spacing is given by Slu = Tm^c/L. The single-mode de- 
scription of Sec. [TV] is valid as long as Suj is much larger 
than all other frequency scales, which poses the restric- 
tion L < 10 cm for mechanical frequencies of lo t /2-k < 
100 MHz. Such a cavity could, e.g., be realized by ter- 
minating an optical fiber with coated end mirrors. For 
the mirrors demonstrated in Ref. [64] (amplitude loss of 
85 ppm) we obtain the naive estimate Ko//27r ~ 10 kHz 
for the decay rate of a cavity with L = 10 cm and assum- 
ing a power-attenuation of 1 dB/km, the intrinsic losses 
are of the same order. 



B. Spin and charge qubits 

In principle, the proposed OMT works for all basic 
types of qubits that can be coherently coupled to a me- 
chanical resonator. In the following, we discuss two 
prominent examples. 

1. Spin qubits 

Prime examples for spin qubits |12j are phosphor 
donors in silicon [13] or the ground state triplet of NV 



centers in diamond [14J . Both systems exhibit excellent 
coherence properties, i.e., T 2 = 1.4 ms for NV-centers 
[55] and T 2 = 10 ms for phosphor donors [13j , with more 
recent experiments even achieving T 2 on the order of sec- 
onds [ST]- Also, coherent control of both systems via 
external fields has been demonstrated. 

The coupling to the mechanical resonator is accom- 
plished by attaching a magnetic tip to the beam at the 
position of an antinode as shown in Fig.|8][b). This tip 
produces a magnetic field at the location of the spin 
which depends on the resonator displacement, and to first 
order in the latter we obtain 

H q = ±B .Z+^(b + tf)a z , (73) 

where a is the vector of Pauli matrices, while Bq is the 
free qubit splitting, into which the constant component 
of the tip's field has been absorbed. The qubit-resonator 
coupling is given by A = g s ^BO-a |V.B| /h, with g s m 2 
being the gyro-magnetic ratio, ^ the Bohr magneton, 
and |V-B| the field gradient produced by the tip at the 
location of the qubit. Provided that the qubit splitting 
can be arranged to be Bq = (ui q ,0, 0) with uj q w w r , 
we may store the qubit in the a^-eigenstates and after a 
relabeling of axes (x f-> z) and subsequent RWA directly 
obtain the coupling of Eq. ([2| , which was the starting 
point for our analysis. 

For qubits whose natural energy splitting u> does not 
match the resonator frequency (i.e. ujq 3> w r ), we can 
apply a classical drive of frequency u>d to bridge the 
gap between the two energy scales. In this case we 
take Bo(i) — (251 cos(wdi), 0, ujq), where f2 is the Rabi- 
frequency. After moving to a frame rotating at u>d we 
drop rapidly oscillating terms and transform to a new 
eigenbasis basis rotated by an angle 9. The latter is given 
by tan0 = fi/A, where A = uiq — cud is the detuning of 
the qubit from the drive. The qubit Hamiltonian then 
becomes 

H q = U ^-a z + ~(b+ 6t) (cos(6)a z - sin(0)a x ) , (74) 

with the effective energy splitting w q = \/A 2 + f2 2 <C c^o- 
Provided that we have A <C u q ~ uj r , we may drop all 
rapidly oscillating terms in a RWA and retain only 

-ffq-rcs = ^ sin(0) ( ( r-6 t + a+b) , (75) 

which is of the form given in Eq.j2} if we redefine 
sin(0)A A. 

We take the magnetic tip coupling the beam to the 
qubit to be fabricated of highly magnetic Co7oFe3o which 
exhibits a magnetization of M — 2.3 T//j,q 35 . For 
a cone-shaped tip of height and radius 100 nm we es- 
timate the field gradient 25 nm below the tip to be 
|VB| = 10 7 T/m. Using a resonator of low frequency 
w r /27r « 5 MHz and thus high zero-point motion a ~ 
1.8 x 10 -13 m, we obtain a coupling of A/27T w 50 kHz. 
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Finally, we estimate the degradation of the toroid's qual- 
ity factor due to the presence of the magnetic tip, which 
is placed a distance d > 2.5 (im away from the rim 
of the toroid as discussed above. We make the con- 
servative assumptions that all light which hits the tip 
gets lost and obtain Q > 27re 2xd V mo de/Aj4 t ip- Here, 
Kiodc = / d 3 r e r (r) \E(r)\ 2 / \E max \ 2 is the mode volume 
of the WGM, Atip the cross-sectional area of the tip and 
X the decay constant of the evanescent cavity field. Us- 
ing the values Af ip ~ 10~ 2 ^m 2 , U mo d c = 90 /Ltm 3 [63] . 
X -1 = 250 nm [5U] we obtain Q > 10 13 which is well 
above the intrinsic Q quoted above. 

Based on these numbers, the results presented in 
Sec. |III C] and Sec. |IV C| show that quantum networking 
operations for spin qubits are feasible, both in the long- 
and short-distance setups. 



object and an estimate for the Q-factor is given by Q > 
e^VmodoAVciectrodoeJ.) [75], where e' r is the imaginary 
part of the relative permittivity of the tip and Ki ec t ro de = 
2.5 xlO -3 /jm 3 its volume. Using ej, < 10 for noble metals 
at optical frequencies [73] we find Q > 10 12 , which again 
exceeds the intrinsic Q. 

As has been discussed in Sec. |III| and |IV[ the fideli- 
ties for quantum networking operations resulting from 
the above numbers are mainly limited by intrinsic qubit- 
dephasing. Improvements of Ti would thus have a big 
positive impact on the performance and also, especially 
superconducting qubits would profit from leaving the 
regime where the evolution of the OM degrees of free- 
dom is adiabatic. 



VI. CONCLUSIONS & OUTLOOK 



2. Charge qubits 

A second promising qubit candidate for the setup de- 
scribed in this work is the superconducting charge qubit. 
It has evolved from the original Cooper pair box qubit 
[66j to more recent designs such as the transmon |671 I68j . 
For the latter, coherence times of T-i = 2 /is have been 
observed [69j and recently, the relaxation times have been 
pushed up to Xi = 200 fis [70]. Control or measurement 
tasks can, for example, be performed in a circuit QED 
architecture |20) . 

Coupling a charge-based qubit to a mechanical res- 
onator is most easily achieved by pushing excess charges 
onto an electrode located on the resonator as depicted in 
Fig.[8^c). The resulting coupling amounts to a position- 
dependent capacitance and has been widely discussed 
and also implemented in experiment |38H41j . However, 
since the intrinsic level splitting ujq of superconducting 
qubits is typically in the GHz range, we need to bridge 
the gap to the MHz resonator frequency u r in order to 
obtain a Jaynes-Cummings coupling. Therefore, we mod- 
ulate the gate voltages as described in Ref. [3U] [71] with 
frequency u>d — ~ <^V and amplitude U r °, which yields 
the coupling of Eq. ([2| with 



Here, Ec is the charging energy of the qubit, C r the 
equilibrium resonator capacitance, and do the equilib- 
rium electrode separation of the coupling capacitance. 
For E c /h w 20 GHz, d « 100 nm, and V r ° a couple of 
volts one easily obtains values in the MHz regime (see 
also Refs. [3J|J HJJ). If it is necessary to place the qubit 
further away from the resonator, the coupling could also 
be mediated by an intermediate normal-conducting wire, 
or additional superconducting circuitry |72j . 

We estimate the optical Q-factor due to the presence 
of the electrode, which we assume to cover one third of 
the resonator and to be thinner than the typical skin- 
depth of ~ lOnm. The field then penetrates the whole 



In conclusion, we have discussed the potential of micro- 
mechanical resonators to mediate interactions between 
dark qubits and light. We have described in detail the 
dynamics of a single OMT and found that the condi- 
tions for low-noise operation overlap with those for OM 
ground-state cooling. Further, we have derived an effec- 
tive description of a multiqubit cascaded network that is 
suitable for long-distance quantum communication, dis- 
cussed the relevant noise-sources in such systems, and 
presented an efficient protocol for state transfer between 
two nodes. In addition, we have shown how OMTs can be 
utilized in an on-chip setting to perform entangling gates 
between qubits that are spaced by less than ~ 10 cm. 
Here, OMTs might provide an alternative to other cou- 
pling schemes [3H31 13Z] for scalable quantum processors. 

We point out that the ingredients needed for the OMT 
are quite generic and may be implemented with a vari- 
ety of solid-state qubits and different OM systems in the 
optical, as well as the microwave domain [73] [75]. In this 
work, we have specifically discussed a potential realiza- 
tion for spin and charge qubits, showing that long- and 
short-distance communication is feasible for present-day 
technology. The gate-schemes presented in this work only 
provide a starting point for optimized and more complex 
protocols. In particular, the long-distance state trans- 
fer protocol could be extended by schemes to correct for 
photon loss (see, e.g., Ref. [75]), while the on-chip gates 
could profit from spin-echo techniques to reduce the im- 
pact of qubit dephasing. In both settings, the limitations 
posed by imperfect qubits could be overcome by depart- 
ing from the adiabtic desciption of the OMT and devis- 
ing fully dynamical gate sequences. This seems especially 
worthwhile for charge qubits, where the qubit-resonator 
coupling A can be made large rather easily. The benefit 
would be shorter gate-sequences at the price of possibly 
stronger requirements on controllability. 

More generally, the OMT can be used to realize com- 
plex qubit-light interactions, which may have applica- 
tions in the implementation of nonlinear optical devices 
on a few-photon level. These are particularly interesting 
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in the context of integrated nanophotonic systems, where 
the proposed device naturally fits in. 
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relations 

c om v = 
vc om = 
rc int r P = o 



(steady state) , 

(Com preserves trace) , 

(vanishing of Tr om {p°™bi} = (b t } {Tee ) , 



where the last conditions follow from the fact that we 
have removed any classical forces from the description 
of the OM systems. These relations allow us to drop 
Com everywhere except for the exponent and since we 
are interested in the dynamics to second order in A we 
may subsequently drop C m t in the exponent leaving us 
with 



Appendix A: Elimination of the optomechanical 
system 

We briefly describe the elimination of the OM degrees 
of freedom to obtain an effective description of the qubit 
dynamics. We denote by p the density matrix of the com- 
plete A^-node system composed of qubit and OM degrees 
of freedom and our goal is to derive an effective ME for 
the reduced qubit density operator p = Tr om {p} to sec- 
ond order in the qubit-resonator coupling A by means of 
projection operator techniques [51]. The way the various 
OM systems interact with each other is left unspecified 
in this Appendix and we lump the whole OM dynamics 
into a Liouvillian £ om equivalent to the linearized OM 
QLEs of the setup under consideration. In an interac- 
tion picture with respect to the free qubit Hamiltonian 
uj q J2j &\ /2 the ME for the complete setup then reads: 



P : 

Cmt(t)p ■ 



C(t)p = 



(Co 



C int (t))p, 
+ b\a; e 



(Al) 
(A2) 



where Cmt describes the qubit-resonator inter acti ons. 
Note that aj ', fo; can be nodal operators (Sec. 



Ill 



and 

App.[B| or normal mode operators (Sec. |IV| and App.[C| 

and that the commutators of the a~ , will not enter 
from now on. We define the projector on the relevant 
part of the density matrix as V p = Tr om {/5} eg) p°™, where 
p°™ is the steady state of the OM system in the absence 
of the qubits, defined by C onl p°™ = 0. We also intro- 
duce Q = 1 - V and note that V 2 = V, Q 2 = Q and 
VQ = QP = 0. To derive an effective equation of motion 



for p we project Eq. ( Al ) on the V- and Q-subspaces and 



formally integrate the equation for the Q-part, where we 
assume for simplicity that at the initial time Qp(t ) = 0. 
The formal expression for the "P-part is then 



Vp = TC(t)Tp 



+ VC(t)Q I dsTexp f &tQC(t)Q 

J to Y.J 8 



QC(s)V P (s) , 



where Texp[. . .] is the time-ordered exponential. Consid- 
erable simplification is brought about by exploiting the 



Vp(t)=VCmt(t) 



At 



Cmt(t-T)V P (t-T) 



(A3) 



Upon expanding the interaction Liouvillian and taking 
the trace implicit in the left-most ^-operation we ob- 
tain terms of the following form (neglecting transients 
by sending to —> — oo) 



p cx A 



dre ± ^ T /(r)£ A £ BM (t-r) 



where /(r) = (A(r)S(0)) fref - Tr om {Ae £ ™ r Sp° s m } is a 
resonator-resonator correlation function evaluated in the 
steady state of the OM system and Ca,b are associated 
qubit Liouvillians. For illustration, we assume that /(t) 
has a dominant contribution f(r) cx e^ %n ~ 1 ^ T and with 
e = i(Q ± uj q ) — 7 integration by parts then yields 

A 2 A 2 f 00 

p oc — CaCbP^) / dr e eT CaCbP^ — t) . 

e 6 Jo 

By iterating this equation we obtain an expansion in 
A/ |e| and for A <C |e| the first, Markovian term domi- 
nates. If this is true for all terms in Eq. (A3) and for 
all frequency components of the correlation functions, 
we may neglect the non-Markovian corrections, which 
is equivalent to replacing p(t — r) — > p(t) in Eq. (A3 1. 
In general, this procedure is thus valid if the OM nor- 
mal modes which couple to the qubits decay much faster 
than A -1 or if they are detuned from the qubits by much 
more than A. The final result can be compactly writ- 
ten as a Markovian ME with coefficients given by one- 
sided Fourier transforms of resonator correlation func- 
tions evaluated at the qubit frequency: 



A 2 r 

fi = --T- 2J l s ij(u q )(afajp - <Tjpa+) 



(A4) 



+T ij (w q )(a i a+p - <r+p<T l ) + h. 



where we have dropped terms rotating at exp[±i2oj g i] 
that contain operators such as crfcr^ based on the as- 
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sumption A <C ui q . The coefficients are given by 

/>oo 

Saiu) = / dr (6 l (r)6 t (0)) frcc e l " T , (A5) 
Jo 

/>00 

T y (w) = / dr (6j(r)6 J (0)) froc e"*" , (A6) 



and for later convenience we also introduce Xij(u) — 
Sij(u)) - T^{u) and Y^lo) = T^u) + T^oj) such that 

p OO 

X ij (u)= dr([bi(T)M(0)}) iiee e iUT , (A7) 



dr(6i(r)6 J (0)) free e-^. (A8) 



It will turn out that interaction- and decay-rates are 
given by X y (w), while Yy (w) occurs in diffusion terms. 



Appendix B: Master equation for cascaded setup 

Here, we derive the effective ME for N qubits cou- 
pled to a chain of cascaded OM systems as discussed 
in Sec. |III| Also, for N = 1 we reproduce the results of 
Sec. [TT] in a different language. To begin with, we rewrite 
the general result ( A4 ) to separate single- and multiqubit 



terms (this rewrite is exact): 



/' = H \ -i-rr 1° -. •/' 



I ^(N t + l)V[a-} t ,+ I ^N l V[at}^ 



where T>\a\\x — 2a^ia i — a 1 a [l — fio^a denotes a Lindblad 
term with jump operator a. Here, we have defined the 
effective frequency shifts Aj = Ao,j + A t h,i with contri- 
butions 

A 2 A 2 
A ,i = — lm{Xu(LO q )} , A t h,i = -—ha.{Tu(u}q)} , 

(B2) 

as well as the effective decay rates I\ and occupation 
numbers iVj, 

I\ = yRe^K)} , N 2 = ^-y«(« t ) . (B3) 

Further, the multiqubit coupling and diffusion rates are 
given by 



_ A 2 A 2 
Ja — —Xij(uj q ) , D.j j 



y — - ^ (B4) 



respectively, and we stress that the cascaded nature of 
the coupling will only become evident below, where we 
evaluate Jy from the physics of the underlying, elimi- 
nated system. 



The OM steady state correlation functions (A5)-(A8) 



determining the effective dynamics are conveniently cal- 



culated from the QLEs ( 30 ) and the input-output relation 



(281 by Fourier or Laplace transformation. To give the 
resulting expressions in a compact form we use the OM 
response matrices A % (i>S) = (M l — iujt)^ 1 and in addition, 
we define the single- node transfer matrices C l (ui) that 
describe how the fiber-field is modified when bypassing 
node i. They are given by C 1 (lu) = (P - 2n l f PA i (uj)P), 
where P = diag[0, 1, 0, 1]. All relevant quantities can 
then be written in terms of the multinode response- 
matrices T IJ (w) defined by (i > j > 0) 

T dj {u) = Ai(u), 

t j h = -2 A /«>'4^ i Hc ,i - 1 (w) • • • c 3+1 ((j)w j (u) . 

In addition, to aid the description of the common fiber 
input, we also introduce (i > 0) 



(B5) 



To determine the correlation functions Xij (ui) we make 
use of the quantum regression theorem |49j with the ini- 
tial condition ([&i(0), &J-(0)])free = $ij an d solve the re- 
sulting equations in the Laplace domain. Since node i is 
only driven by nodes j < i, we obtain 



XiAw) 



for i < j 

7ft (w) fori>j 



(B6) 



which, via Eq. (B4), turns Eq. (Bl I into a cascaded ME 



In order to specify the noise terms in the ME we 



(Bl) 

solve Eq. (j30h in the frequency domain, with the Fourier 



transformation of some quantity f(t) defined as f(u>) — 
fdte tuJt f(t) (which induces [/(w)]t = f(-u>)). As 

a result, we get v'H = 7 i V) I M-Ej<i T ij (u)B. j (u) , 
where I(w) = I 1 ^) = (0, f ia (u), 0, f in (u)) T is the fiber 
input of the first node. Due to the 5-correlated nature of 
the noise sources we can express the second moments of 
the OM operators as follows = 1 . . . N label nodes 
and k, I = 1 . . . 4 vector components): 



(<VWM)free = C%(u)S(u' + Ul) , 



(B7) 



Here, we have introduced the matrices C lJ (w) , which di- 
rectly give the desired correlators for the ME via Yij [ui) — 
C^(w). They can be expressed in terms of the transfer- 
matrices and noise statistics according to 

C^{oj) = {T i0 (lu))*I°(T i0 (uj)) t 

m'm(i.j) 

+ £ (r m H)*r"(r j "(^)) T , (B8) 
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where the matrices r l characterize the noise correlations 
of node i by means of (B% (t)Rj(t')) = r l kl 5{t - t'), while 

1° describes the fiber input, (4 (*)!;(*')) = ^A* ~ *')■ 
From the statistics quoted in the main text it follows 
that r* = diag[ 7 ^7V^,0, 7 ^(7V; n + l),2«o] and 1° = 
diag[0, 0,0,1]. 

Finally, we have to evaluate the thermal frequency shift 
Ath.i, which we do exemplarily for the first node. Us- 
ing the quantum regression theorem we obtain the rele- 
vant correlation function Tn(u>) = ■^■ik( ul )( v k v l)free> 
where the steady state moments are given by the fre- 
quency integral (v k 'vl)f TCC — J |^C^(cj), which can be 
evaluated exactly [52] . If the qubit is on resonance with 
a normal mode of the OMT we obtain for A r = u r 



A 



7i 



,.N„ 



G 2 A 



K 

2wF 



(B9) 



where we have neglected higher orders in K/u> r , G juj r as 
well as non-RWA corrections to the thermal part. From 
the low noise conditions identified in Sec. Ill El it follows 
that | A t h | "C |Ao|. If, however, the qubit is off-resonant 
with the OM system as for G — > on the trajectory of 
Fig.[2jc), we find that the leading orders are still given 
by the above expression. In conclusion, we thus have 
| A t h | ;$ |Ao| for most cases of interest and finally note 
that A t h can be compensated to the same amount as the 
unavoidable shift Ao by tuning of the bare qubit frequen- 
cies. 



Appendix C: Master equation for setup with fiber 
cavity 

Here, we derive the effective qubit dynamics for the 
case of the on-chip setup discussed in Sec. |lV[ Note that 
the results given here apply to the case of identical nodes. 



Network with N identical nodes 



According to the second term in Eq. ( 60 ) , only the 
" center-of-mass" mode couples to the fiber cavity and 
we therefore introduce normal modes for the cavities ac- 
cording to 5j = VijCj with V±j = 1/ V~N, and similarly for 
mechanical resonators and qubits. The remaining rows 
of V are constructed from the constraint V T V — 1, such 



that the remaining terms in Eqs. (591,(60) are left invari- 
ant. As a result, we can write the QLEs describing the 
OM variables as N decoupled sets 



y(t) = -M i v i (t)-~R i (t), i = l. 

(b 1 ,c 1 ,d ,b\,c\,dl) T 



,N . 



(CI) 



,%i>2 



and V 

(bi, 6j, b\, c\) T . Note that the systems i > 2 
are all identical, which will simplify the final re- 
sults. The noise vectors R l (i) can be constructed 
easily from the original QLEs (57) and the or- 



they are mutually uncorrelated. We have R 1 (i) = 

(V^mfl) V2Ko/o,l, V / 2Ko//o, VTmih V^ofo,l> V^QffoV 

and R^ 2 (t) = (V7^^,y2^/ 0il ,77^e7,\/2^/o, l ) T - 
The statistics of these operators are given by 
((Rl(t)Y Ri(t')) = ? l kl 5(t - t') with diagonal correlation 
matrices f 1 = diag[7 m iV ro , 0, 0,j m (N m + 1),2k ,2k 0/ ] 
and f 1 ^ 2 = diag[7 m A^ m , 0, j m {N m + 1), 2k ]. Finally, the 



drift matrices of the systems (CI ) are M l ~ 2 — M, where 
M is defined in Eq. ( 15 ) with k — > kq and 



M 1 = i 



( 


L0 r 


G* 








CG 


\ 




G 


A c 


K* 


CG 













K 


A c0 
















-CG* 





— U3 T 


-G 















-G* 


-A c 


—K 


V 














-K* 


-A c0 / 



D. 



where D = diag[7 TO /2, kq, K /,7m/ 2 , «o, «o/]- Here, K = 
i/Nh is the coupling of the center-of-mass mode to the 
bus cavity and the full OM coupling is again obtained for 
C = 1, while in RWA the off-diagonal blocks are dropped 

(C = o). 

We proceed to eliminate the OM degrees of freedom 
by applying the results of App.[S] based on the inter- 
action H lnt = ^^vivtbv + h.c.)/2 between qubits and 
resonators. The N sets of normal modes described by 
Eq. (CI I are mutually uncorrelated, i.e., ( 6 j (t ) St (0) ) f ree oc 



5ij, and we may hence eliminate them one by one. There- 
fore, we drop the superfluous second index on the cor- 
relation functions (A5)-(A8) and finally note that the 
correlators for sets i = 2 . . . N are all identical, that is, 
Xi> 2 {<jj) = X 2 (uj), etc. The effective ME (fA4| may then 



be written in terms of collective and local contributions: 



(i=-i[j N S + S - i— [S z ,n] 



(C2) 



^{r coU (iv col i + i)v[s-]fx + r coU N con v[s+]fi} 
i N 

o 12 { r ioc(Moc + l)V[<J-}ti + T loc N loc V[a+]^ 



Here, we have introduced the collective operators S = 
and S z — X^j^zj an d the interaction and decay 
rates are given by 

J w = ^Im{A > 1 -A > 2 }, A=^Im{X 2 }, (C3) 

Tcoii = ^ Re {^i - ^2} , r loc = yRe{X 2 } , (C4) 

where we neglected a thermal contribution to A. The 
collective and local occupation numbers are determined 
by 



thogonality of the transformation V ensures that 



A2_ 
AN 



(y x - %) , r loc 7v loc 



A 2 - 

-rYi 



(C5) 
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and all correlation functions are understood to be eval- 
uated at the qubit frequency uj q . They can be eval- 
uated as for the cascaded setting and with A l {w) = 
(M* - iujty 1 we obtain Xi(u) = Aj x (o;) and Yi(u>) = 
[l i *(w)f i i iT (u;)]ii. 



where A' = A + J with J = J2 and A given above, and 



T v = ~Re{X v } , VN V = ^-Y v . (C7) 



Two node network 



For N = 2 identical nodes the sets v 1 and v 2 corre- 
spond to the symmetric and anti-symmetric sets intro- 
duced in the main text as v s and v a , respectively (cf. 
Eq. (61 1). In this case, the ME (C2) simplifies since the 
Lindblad terms are diagonal in the symmetric and anti- 
symmetric qubit operators a~ a = (cr£~ ± o^~)/v2: 

A' 

fi = -i[j{a^a^ +crfo-^),/i] - iy ^ [<?l,(j] (C6) 

i=l,2 



Finally, we estimat e the fi delity for generating a Bell state 

We add a term j^- ^[^z] to 



as discussed in Sec 



IV C 



Eq. ( C6 1 in order to account for intrinsic qubit dcphasing 



and obtain to first order in the Lindblad terms: 
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